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SUMMARY 


Research  Goals 

The  elastohydrodynamic  (EHD)  regime  of  lubrication  is  responsible  for  the  generation  of  a 
protective  film  and  the  transfer  of  shear  across  concentrated  contacts  in  many  machine  elements 
including  bearings,  cams  and  gears.  It  is  now  understood  that  small  scale  EHD  plays  a  major  role 
in  the  boundary  lubrication  of  rough  surfaces.  Many  high-pressure  metal  working  operations 
share  the  pressure  and  kinematics  of  elastohydrodynamics.  A  complete  solution  of  the  EHD 
problem,  however,  requires  a  thorough  knowledge  of  the  lubricant  constitutive  behavior  and  the 
attending  properties  as  functions  of  temperature  and  pressure.  The  Newtonian  assumption  alone 
is  inadequate. 

It  has  been  observed  in  boundary  lubrication  and  in  some  aspects  of  elastohydrodynamic 
lubrication  that  friction  is  nearly  Coulombic  in  nature  -  the  fiiction  coefficient  is  only  weakly 
dependent  upon  load  and  sliding  velocity.  In  some  instances  the  fnction  coefficient  may  be  so 
similar  in  the  boundary  and  EHD  regimes  that  fiiction  alone  does  not  clearly  discriminate  the 
transition  from  one  to  the  other.  These  attributes  of  liquid  lubrication  would  seem  enigmatic. 

It  has  been  the  goal  of  this  program  to  elucidate  the  deformation  mechanisms  responsible  for  the 
unexpected  behavior  of  liquid  lubricants  under  the  uniquely  severe  conditions  of  EHD  and  provide 
numerical  analysts  and  designers  with  a  rheological  model  for  prediction. 

Significant  Results 

Work  supported  by  ONR  and  performed  in  this  laboratory  has  resulted  in  major  advances 
in  understanding  the  deformation  mechanisms  of  lubricants  at  high  pressure.  Among  these  are: 

1)  Rate  independent  and  perhaps  even  non-linear  shear  rheological  response  for  liquids  at 
high-pressure  and  high-stress  is  the  result  of  intermittent  slip  along  inclined  planes  which 


are  consistent  with  the  Mohr-Coulomb  theory  if  a  first  normal  stress  difference  is 
assumed. 

2)  Slip  within  the  lubricant  fiJm  is  associated  with  dilatation  which  results  in  the  pressure 
dependence  of  the  shear  strength  of  these  films. 

3)  Deformation  within  a  lubricant  film  may  localize  by  a  thermal  mechanism  as  well,  if  the 
Brinkman  number  is  suflSciently  large  and  the  flow  is  stress-controlled. 

4)  The  pressure  capability  of  high-pressure  high-stress  rheometers  has  been  improved  so 
that  accurate  calculations  of  concentrated  contact  traction  can  be  obtained  fi'om  primary 
material  property  measurements. 

Plans  for  Next  Year 

We  have  introduced  the  Mohr-Coulomb  failure  criterion  as  a  method  for  predicting  slip. 
Mohr-Coulomb  defines  two  possible  orientations  for  stress  induced  shear  bands.  Both  are 
experimentally  observed  in  our  high-pressure  flow  visualization  cell  and  the  measure  of  the 
included  angle  between  the  types  of  bands  is  consistent  with  theory.  A  first  normal  stress 
difference  must  be  introduced  to  account  for  the  orientation  of  the  shear  bands  with  respect  to  the 
principal  shear  directions. 

The  first  normal  stress  difference  in  uniaxial  shear  was  the  subject  of  conjecture  and 
controversy  early  in  the  study  of  high-pressure  rheology  of  lubricants;  however  there  has  been  no 
attempt  to  measure  this  effect  under  high-pressure.  An  existing  high-pressure  rheometer  will  be 
modified  so  that  the  normal  stress  may  be  measured.  Shear  stress  and  the  normal  stress  difference 
will  be  reported  as  functions  of  the  shear  rate  (i.e.,  the  \nscometric  functions). 


4 


The  pressure  capability  of  our  flow-visualization  system  has  been  increased  to  1  GPa  in  a 
squeeze  film  configuration.  We  will  extend  this  new  pressure  capability  to  the  unaxial  shear 
configuration  for  which  the  interpretation  of  results  is  unambiguous.  We  can  then  seek  a 
correlation  between  the  nucleation  of  shear  bands  and  the  viscometric  functions. 

The  viscometric  functions  for  shear  and  normal  stress  differences  will  be  incorporated  into  a 
model  for  shear  band  prediction.  The  consequences  of  a  successful  program  are  far-reaching. 
The  fiiction  coefficient  in  concentrated  contact  has  been  treated  as  a  disposable  parameter  in 
numerical  and  analytical  studies.  We  will  show  that  the  fiiction  coefficient  is  instead  related  to 
material  properties  of  the  lubricant. 
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1.  BACKGROUND 


The  elastohydrodynamic  (EHD)  regime  of  lubrication  is  responsible  for  the  generation  of  a 
protective  film  and  the  transfer  of  shear  across  concentrated  contacts  in  many  machine  elements. 
It  is  now  understood  that  small  scale  EHD  plays  a  major  role  in  the  boundary  lubrication  of  rough 
surfaces.  Many  high-pressure  metal  working  operations  share  the  pressure  and  kinematics  of 
elastohydrodynamics.  A  complete  solution  of  the  EHD  problem,  however,  requires  a  thorough 
knowledge  of  the  lubricant  constitutive  behavior  and  the  attending  properties  as  functions  of 
temperature  and  pressure.  The  Newtonian  assumption  alone  is  inadequate. 

Early  investigations  of  high-pressure  lubricant  rheology  addressed  the  small  strain 
response.  For  small  strains,  viscous  heating  which  plagues  large  strain  measurements  may  be 
ignored.  These  measurements  were  both  optical  and  acoustical,  e.g.  references  [1]  and  [2] 
respectively,  and  provided  linear  viscoelastic  properties.  These  techniques  are  experimentally 
convenient  in  requiring  only  optical  or  electrical  paths  through  the  high-pressure  chamber.  The 
linear  viscoelastic  response  is,  however,  unable  to  account  for  any  but  the  small  slip  traction  in 
EHD. 

Attention  then  turned  to  large  strain  investigations  which  could  reveal  non-linear  behavior 
and  successfully  predict  EHD  traction  [3].  The  progress  in  this  area  has  recently  been  significant 
in  that  the  mechanism  of  non-linearity  (at  least  for  low  molecular  weight  liquids)  has  been  shown 
to  be  a  mechanically  induced  shear  localization  or  shear  bands  [4].  The  constitutive  behavior 
which  in  a  rheological  measurement  by  definition  excludes  the  shear  bands  is  apparently  linear 
viscoelastic.  Attention  must  then  eventually  return  to  the  small  strain  linear  viscoelastic  behavior. 
For  it  is  the  constitutive  response  which  provides  the  stress  field  which  initiates  the  localization. 
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That  is  to  say,  one  must  understand  the  constitutive  response  of  the  homogeneously  deforming 
matrix  to  predict  the  onset  and  orientation  of  shear  bands. 

Fundamental  to  the  generation  of  a  film  in  concentrated  contact  is  the  piezoviscous 
property  of  liquid  lubricants  -  whereby  even  simple  low  molecular  weight  liquids  at  ambient 
temperature  attain  very  great  viscosities  under  pressure.  Shear  rheological  investigations  of  non- 
Newtonian  response  of  lubricants  [5,6]  have  resorted  to  low  temperatures  to  achieve  the  same 
level  of  viscosity.  It  is  clear  now,  that  high-pressure  is  essential  to  an  accurate  simulation. 
Liquids  become  considerably  "stronger"  under  pressure.  Peculiar  to  these  atmospheric  pressure 
studies  was  the  observation  of  fi-acture.  Eastwood  and  Harrison  [5]  could  observe  the  liquid  in 
shear  and  reported  cracking. 

Early  investigations  of  non-Newtonian  lubricant  flow  often  studied  polymer  solutions. 
The  blend  of  mineral  base  oil  and  polymeric  viscosity  index  improver  is  representative  of  multi¬ 
grade  motor  oils.  The  capillary  rheometer  was  useful  for  this  purpose  as  it  is  very  simply 
pressurized  and  although  the  shear  stress  caimot  be  large  for  this  instrument,  non-Newtonian  flow 
occurs  at  relatively  low  stress  for  polymer  solutions.  The  high-pressure  capillary  viscometer  was 
developed  to  a  high  level  by  Jakobsen  and  Winer  [7]  who  reported  measurements  to  shear  rates  of 
10  VI  ,  pressures  to  600  MPa  and  time  of  shear  as  short  as  4  ps.  They  reported  Newtonian  flow 
for  liquid  lubricant  base  stocks  to  a  shear  stress  of  5  MPa. 

A  flow  chart  for  a  mineral  oil  /  PAMA  blend  is  reproduced  fi'om  [8]  in  Figure  1  for 
pressures  and  temperatures  indicated.  The  logarithm  of  effective  viscosity,  p*  =  t/y,  is  plotted 
against  the  logarithm  of  shear  stress,  x.  The  diagonal  lines  represent  rate  of  shear,  y .  The  typical 
transition  fi'om  the  first  to  the  second  Newtonian  regimes  are  evident.  The  liquid  flows  with 
constant  effective  viscosity  provided  that  the  stress  is  below  or  above  the  transition  region.  The 
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presentation  in  Figure  1  is  particularly  effective  in  that  the  non-Newtonian  transition  is  shown  to 
be  defined  by  a  more-or-less  constant  critical  stress  rather  than  a  critical  rate.  That  is,  individual 
flow  curves  may  be  superimposed  by  shifting  along  a  linear  path  which  to  an  approximation  is  a 
constant  shear  stress  path.  This  observation  (which  has  not  always  been  well  recognized  in  the 
Tribology  community)  is  a  statement  of  time-temperature  (-pressure)  superposition.  Following 
reference  [9],  to  obtain  the  effective  viscosity  at  a  temperature,  T,  fi'om  the  flow  curve  measured 
at  a  reference  temperature,  Tr,  set 

|i.(T.Y)  =  ;^l*.(T,.a,Y)  (1) 


where  p  is  the  limiting  low  shear  viscosity  and  aj  is  the  shift  factor.  A  good  approximation  for 
ax  is  [9,10] 


"  ~  KTr) 


(2) 


The  superposition  technique  was  extended  to  include  pressure  by  Phillippoff  [1 1]. 
Equations  (1)  and  (2)  yield  the  result 

P.  =mP 


(3) 


where  P  is  the  temperature- viscosity  coefficient,  den\i/dT,  and  m  is  the  rate  sensitivity  coefficient, 
diti'c  I  deny ,  and  pg  is  the  effective  temperature-viscosity  coefficient  for  a  rate-controlled  process 
with  non-Newtonian  response.  Equation  (3)  was  also  found  experimentally  by  Winter  [12]  and  is 
very  useful  in  assessing  viscous  heating  effects  for  high  dissipation  rheometry.  It  should  be  clear 
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Log  Shear  Stress/Pa 


Flow  Chart  for  Blend  of  Parafllnic  Mineral  Oil  +  4% 

Polyalkylmethacrylate  from  [8]. 


from  Figure  1  and  the  above  that  a  rheological  experiment  to  obtain  constitutive  information  for 
lubricants  must  simulate  the  shear  stress  of  the  lubrication  problem  rather  than  the  rate  when  the 
dissipation  is  limited  by  experimental  technique.  In  addition,  an  accurate  correlation  of  the 
limiting  low  shear  viscosity,  p,  with  temperature  and  pressure  is  indispensable  to  an  investigation 
of  constitutive  behavior  of  liquids  under  pressure. 

2.  VISCOSITY 

The  measurement  of  limiting  low  shear  viscosity  (hereafter,  simply  viscosity)  under 
pressure  has  become  relatively  routine.  Some  commercial  viscometers  can  reach  300  MPa. 
Laboratory  high-pressure  viscometers  can  operate  to  1  GPa  and  above  and  are  most  often  of  the 
falling  body  type  which  is  described  in  [13].  Falling  body  viscometers  apply  a  shear  stress 
typically  less  than  100  Pa  so  that  the  liquid  response  should  be  NeAvtonian;  although  the  rate  of 
shear  may  exceed  lO^s'l  for  low  viscosity.  Because  of  the  low  applied  shear  stress,  the  sinker  fall 
times  become  inconveniently  large  for  viscosity  greater  than  10^  Pa-s. 

A  new  variation  of  the  falling  body  viscometer  was  developed  in  this  laboratory  and  was 
utilized  to  generate  some  of  the  results  for  viscosity  greater  than  10^  Pa-s  in  Figure  2.  The  sinker 
would  conventionally  be  made  to  fall  under  the  force  of  gravity  toward  the  bottom  of  a  pressure 
chamber.  The  time  of  fall  yields  viscosity.  If  one  places  a  spring  between  the  sinker  and  the 
closed  bottom  end  of  the  pressure  chamber,  any  increase  in  pressure  which  reduces  the  specific 
volume  of  the  liquid,  will  push  the  sinker  toward  the  closed  end  of  the  chamber,  thereby 
compressing  the  spring.  The  spring  will  eventually  return  the  sinker  to  its  equilibrium  rest 
position  while  displacing  liquid  through  the  clearance  (typically  100  pm)  between  the  sinker  and 
the  chamber.  The  system  relaxes  to  equilibrium  with  exponentially  decaying  motion  as  a  spring- 
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the  chamber.  The  system  relaxes  to  equilibrium  with  exponentially  decaying  motion  as  a  spiing- 
dashpot.  The  time  constant  for  this  motion  will  yield  a  measure  of  viscosity  if  the  initial  sinker 
position  is  defined  at  a  sufGciently  late  time  so  that  the  heat  of  compression  has  been  conducted 
fi'om  the  liquid  into  the  chamber  walls. 

The  maximum  viscosity  which  can  be  measured  by  the  falling  body  technique  is  limited  by 
the  stability  of  position  sensing  instruments  over  long  fall  times.  The  spring-dashpot  technique 
increases  the  maximum  obtainable  viscosity  for  a  given  observation  time  by  increasing  the  applied 
shear  stress  to  about  600  Pa.  An  additional  benefit  of  the  increased  stress  is  that  weak  wax 
structures  (which  sometimes  form  in  mineral  oil  at  high  pressures  and  arrest  the  fall  of  a  sinker) 
are  broken. 

The  viscosity  under  pressure  for  two  perfluorinated  polyethers  is  shown  graphically  in 
Figure  2.  The  143AD  and  143AC  diflfer  only  in  molecular  weight  -  the  former  being  higher  in 
molecular  weight.  The  slope  of  an  isotherm  is  related  to  the  local  pressure  viscosity  coefficient, 
a  =  dtn\i  /  dp.  Clearly  the  slope  and  a  decrease  with  pressure  initially.  Then  at  some  inflection 
they  begin  to  increase  with  pressure.  Bridgman  [14]  suggested  that  for  pressure  greater  than  1.2 
GPa,  a  increases  with  pressure  for  all  liquids.  This  was  supported  by  measurements  to  6  GPa  on 
some  organic  liquids  by  Barnett  and  Basco  [15],  and  more  recently  at  lower  pressures  [16], 

Several  useful  pressure-viscosity  correlations  are  avdlable  for  low  pressures  characteristic 
of  the  EHD  inlet  zone.  These  generally  treat  the  pressure-viscosity  coefficient  as  monotonically 
decreasing  with  pressure  and  for  most  liquid  lubricants  are  suitable  to  0.3  GPa  provided  the 
viscosity  does  not  exceed  about  10^  Pa  s.  The  use  of  these  correlations  for  extrapolation  to  EHD 
Hertzian  zone  pressures  can,  however,  yield  misleading  trends.  For  example,  Hsiao,  et  al  [17] 
fitted  the  parameters  of  the  isothermal  Roelands  equation  to  the  viscosity  of  143AC  up  to  0.3 
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Figure  2.  Pressure- Viscosity  plot  for  two  perfluorinated  polyalkylethers.  Solid 
curves  are  Free  Volume  Model  [18].  Broken  curve  is  extrapolation 
from  [17].  Points  are  Measured  Viscosities. 
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GPa  and  used  this  correlation  to  extrapolate  to  1.0  GPa  pressure.  The  broken  curve  in  Figure  2  is 
their  result.  Clearly,  at  only  0.6  GPa  the  measured  viscosity  diverges  from  the  extrapolation. 
Although  this  deficiency  has  been  known  for  some  time  and  was  acknowledged  by  Roelands  [3], 
the  Roelands  equation  continues  in  use  for  numerical  simulation  of  the  lubricant  in  the  Hertzian 
zone. 

The  viscosity  of  the  mineral  oil,  HVI  650,  was  obtained  using  a  falling  body  viscometer 
and  a  Couette  viscometer  which  will  be  described  later.  The  curves  in  Figure  3  and  the  solid 
curves  in  Figure  2  represent  the  Free  Volume  Model  of  Yasutomi  [18],  The  glass  transition  is 
taken  to  be  an  isoviscous  state  with  viscosity  of  Pg  and  with  the  glass  transition  temperature 
related  to  pressure  as 

T,=T,„  +  A/n(l  +  A2p)  (4) 


The  relative  free  volume  expansivity  is  also  taken  to  be  pressure  dependent  as 

F  =  l-B/n(l+B2p) 

The  viscosity  can  be  calculated  from  the  WLF  equation. 


p  =  Pg  exp 


-2.3C,(T-Tg)F 


(5) 


(6) 


The  parameters  for  equations  (4-6)  were  regressed  from  viscosity  data  and  are  tabulated  for  some 
of  the  lubricants  used  in  this  paper  in  Table  1  where  a  value  of  Pg  was  assumed.  Use  of  this  model 
provides  some  contact  between  the  piezoviscous  effect  and  state  properties. 
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Figure  3.  Pressure- Viscosity  plot  for  HVI  650.  Curves  are  for  Free  Volume 

Model. 
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Table  1. 

Free  Volume  Parameters 


Liauid 

Uo/Pas 

A/C 

A,/GPa-‘ 

B,/GPa^ 

Cl 

C7/“C 

143AC 

lO’ 

-64.93 

185.8 

1.628 

0.1806 

37.55 

10.802 

50.69 

143AD 

lO’ 

-60.91 

183.8 

1.805 

0.1277 

75.42 

10.592 

49.44 

HVI650 

lo’ 

-30.4 

309.1 

0.3064 

0.2186 

29.99 

10.264 

27.04 

Santotrac  40 

lo’ 

-79.62 

91.86 

3.805 

0.3865 

5.138 

12.381 

52.14 

SHFlOOl 

10’ 

-58.2 

325.6 

0.2948 

0.1536 

30.79 

11.120 

57.88 

Table  2. 

Comparison  of  Relative  Thermal  Expansivity, 

e(p)/8(o),  with  that  Calculated  from 

Free  Volume  Model  for  143AC 

from  [19] 

Pressure/MPa 

e(p>/"C‘ 

efoVEfol 

F/f2-Fl* 

0 

9.7  X  10^ 

1.00 

1.00 

34.5 

6.9 

0.72 

0.74 

69 

5.7 

0.59 

0.62 

103.5 

5.0 

0.52 

0.55 

138 

4.5 

0.46 

0.50 

172 

4.14 

0.43 

0.47 

*F  from  eqn  (5) 
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Recent  numerical  simulations  of  thermal  compressible  EHD  have  incorporated  state 

equations  which  treat  the  thermal  expansivity,  e,  as  a  constant  with  a  value  equal  to  the 

atmospheric  pressure  value.  This  must  be  suspect  as  it  can  be  seen  that  the  thermal  expansivity 

decreases  to  about  one-third  of  the  atmospheric  value  at  1  GPa  for  most  liquid  lubricants. 

Yasutomi  [18]  showed  that  the  relationship 

e(p)  ^  F(P) 
e(o)  2-F(p) 

holds  for  liquids  of  various  chemical  composition  to  a  pressure  of  at  least  0.6  GPa.  Wilson  [19] 
measured  the  density  as  a  function  of  temperature  and  pressure  for  143AC  using  a  bellows  type  of 
piezometer.  The  thennal  expansivities  were  calculated  from  both  [19]  and  equation  (7)  and  are 
compared  in  Table  2.  The  agreement  should  be  considered  good  in  light  of  the  difference  in 
experimental  technique. 

3,  ELASTIC  SHEAR  MODULUS 

A  description  of  the  transient  small  strain  response  of  liquid  lubricants  at  high-pressure 
requires  as  a  material  property  the  limiting  high-frequency  elastic  shear  modulus  (hereafter  simply 
shear  modulus).  It  was  shown  (Ref  [3])  that  a  Maxwell  element  adequately  represented  the 
transient  uniaxial  shear  stress-strain  response  in  a  rheometer  for  a  mineral  oil  at  a  pressure  of  780 
MPa.  For  small  strains,  only  the  viscosity  and  shear  modulus,  G,  are  required. 

There  have  been  few  recent  experimental  measurements  of  the  shear  modulus  of  liquid 
lubricants  under  pressure.  There  does  exist,  however,  a  body  of  literature  which  shows  consistent 
results  for  the  variation  of  G  with  pressure  near  room  temperature.  These  data  are  collected  in 
Figure  4  where  the  results  of  three  different  techniques  -  light  scattering,  acoustic,  and  direct 


16 


Diester-Barlow 


17 


Figure  4.  The  Elastic  Shear  Modulus  of  Liquid  Lubricants  at  Room 
Temperature  by:  light  scattering  by  Bezot  [20],  acoustics  by  Barlow 
[2]  and  Hutton  [21],  mechanical  by  Bair  and  Winer  [3],  impact  bv 
Ramesh  [22]. 


mechanical  measurement  -  fit  nicely  within  the  envelope  shown  by  the  broken  curves.  References 
are  given  in  the  figure.  For  high  pressures  it  would  appear  that  we  can  approximate  G  »  3p.  Only 
the  pressure-shear  plate  impact  measurements  are  inconsistent.  In  the  impact  work,  thermal 
softening  was  neglected  and  the  observed  non-linear  shear  response  was  attributed  to  constitutive 
behavior  in  the  form  of  the  Eyring  model.  However,  the  authors  estimate  the  temperature  rise 
due  to  shear  heating  to  be  as  high  as  78°C  [22].  In  as  much  as  the  liquid  (5P4E)  studied  is  one  of 
the  most  temperature-sensitive  known,  it  is  more  likely  that  the  non-linearity  can  be  attributed  to 
shear  heating. 


4.  COMPRESSIBILITY 

Lubricant  compressibility  is  a  second  order  effect  for  the  film  thickness  in  concentrated 
contact.  However,  dilatation  apparently  plays  a  key  role  in  the  initiation  of  mechanical  shear 
bands  [23]  which  limit  the  shear  strength  of  liquid  films  as  evidenced  by  the  pressure  dependence 
of  the  shear  strength  [24].  Numerical  simulations  of  EHD  usually  incorporate  a  liquid  equation  of 
state  which  assumes  a  maximum  limiting  density.  This  concept  is  not  supported  by  experimental 
results. 

The  compressibility  of  liquids  is  usually  determined  by  one  of  two  techniques.  The 
displacement  of  a  piston  under  a  known  force  jnelds  a  volume  change  if  the  deformation  of  a 
piston/cylinder  device  can  be  accurately  accounted  for.  Alternatively,  a  piezometer  may  be 
immersed  in  a  fluid  medium  for  which  the  pressure  is  controlled.  Here,  the  piezometer  signal 
must  be  passed  through  the  vessel  which  confines  the  fluid  medium. 

A  unique  instrument  was  designed  and  constructed  in  this  laboratory  and  is  pictured  in 
cut-away  view  in  Figure  5.  The  Differential  Scanning  Dilatometer  (DSD)  is  essentially  two 
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piston/cylinder  devices  placed  in  series  so  that  the  piston  forces  are  identical.  If  the  cylinder  bores 
and  seals  are  identical,  the  pressures  in  the  cylinders  are  equal.  Each  cylinder  communicates  with 
a  separate  compression  chamber  in  a  temperature  controlled  block  through  small  bore  tubing  and 
fittings.  Temperature  of  the  block  is  measured  by  a  thermocouple  centered  between  the  chambers 
and  the  temperature  is  controlled  by  passing  either  liquid  N2  boil  off  or  heated  air  through  a  series 
of  holes  (not  shown  in  Figure  5)  in  the  block.  The  intensifier  vessel,  tubing  and  compression 
chamber  are  filled  with  a  liquid  high  pressure  medium.  The  high  pressure  medium  is  the  diester, 
bis  (2-ethyl  hexyl)  sebecate,  which  is  widely  used  for  high  pressure  work.  Each  chamber  accepts 
a  thin  plastic  capsule  -  one  capsule  containing  a  tungsten  slug  (in  the  reference  chamber)  and  one 
capsule  containing  an  equal  volume  of  the  liquid  sample  (in  the  sample  chamber). 

Two  potentiometers  are  utilized  as  piston  displacement  transducers.  (They  cannot  be  seen 
in  the  Figure.)  The  transducer  outputs  are  electrically  subtracted  so  that  the  signal  represents  xr  - 
xg.  Now,  as  the  system  pressure  is  varied,  the  relative  piston  displacement  xg  -  xr  is  a  measure 
of  the  dilatation  of  the  sample  relative  to  that  of  tungsten.  The  deformation  of  the  vessel  and 
structure  is  removed  fi'om  the  measurement  to  the  extent  that  these  deformations  affect  both 
piston  displacements  equally.  Tungsten  was  chosen  as  a  reference  material  for  its  very  large  bulk 
modulus  and  low  thermal  expansivity. 

Here  we  are  concerned  with  the  isothermal  compressibility.  Indeed,  it  would  probably  not 
be  possible  to  make  measurements  with  sufficient  speed  to  approach  the  adiabatic  idealization. 
To  ensure  that  the  compression  of  the  sample  can  be  regarded  as  isothermal,  a  transient  heat 


19 


INTENSIFIER  VESSEL 


LOW  PRESSURE 
MEDIUM 


PISTON 


HIGH  PRESSURE  MEDIUM 


PISTON 


TEMPERATURE 
CONTROLLED  BLOCK 


REFERENCE  CHAMBER 


SAMPLE  CHAMBER 


Figure  5.  DifTerential  Scanning  Dilatometer  in  Cross-Section. 
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transfer  numerical  analysis  was  performed  to  identify  the  maximum  pressurizing  rate  (which  is  140 
MPa/minute)  for  which  the  heat  of  compression  can  be  effectively  removed  by  conduction. 
Experimentally,  one  could  simply  remain  at  elevated  pressure  for  a  very  long  time;  however  to  do 
so  increases  the  risk  of  leaks  becoming  significant. 

The  instrument  described  above  provides  a  refined  technique  for  observing  liquid 
compressibility  in  a  pressure  range  which  brackets  the  glass  transition  pressure,  pg,  while 
maintaining  a  nearly  hydrostatic  state  of  stress.  For  example,  with  a  piston/cylinder  type  of 
dilatometer,  in  addition  to  reducing  the  volume  of  the  sample,  one  is  also  performing  a  uniaxial 
compression  of  the  cylindrical  sample  and  therefore  introducing  deviatoric  stress.  Also,  shear 
stresses  are  introduced  in  the  interface  between  the  liquid  and  the  cylinder  bore.  Given  sufiBcient 
time,  the  liquid  will  relax  to  a  purely  hydrostatic  condition;  however,  if  the  pressure  is  sufficiently 
greater  than  pg  the  relaxation  time  will  be  inconveniently  long.  The  material  response  using 
previous  techniques  should  be  overly  stiff  for  pressure  greater  than  pg.  The  present  device 
compresses  the  liquid  sample  within  an  easily  deformed  capsule  which  is  in  turn  surrounded  by  a 
low-viscosity  liquid. 

The  accuracy  of  the  present  device  was  tested  by  measuring  the  relative  density, 
p(/3,r)/p(o,r),  of  water  and  the  results  are  compared  with  others  from  the  literature  in  Table  3. 

The  agreement  is  excellent.  New  results  for  the  traction  fluid,  Santotrac  40,  are  presented  in 
Table  4.  The  relative  density  of  5P4E  (long  the  subject  of  high-pressure  studies)  versus  pressure 
is  graphically  depicted  in  Figure  6.  The  glass-transition  is  evident  as  the  "knee”  in  the  curve  and 
agrees  quite  well  with  the  empirical  formula  for  Tg(p)  given  by  Yasutomi  [18].  However,  the 
bulk  modulus  of  the  glass  is  here  considerably  lower  than  that  measured  in  a  piston  cylinder 


21 


device.  The  bulk  modulus  of  glassy  5P4E  is  in  Figure  6  about  7  to  8  GPa  versus  more  than  20 
GPa  for  piston/cylinder  techniques. 

Hamrock,  et  al  [27]  utilized  a  piston/cylinder  type  of  dilatometer  to  report  the  change  of 
density  and  glass  transition  pressure  of  liquid  lubricants  to  2.2  GPa.  Ramesh  [28]  utilized  a 
compression  Kolsky  bar  and  plate  impact  device  to  find  the  compressibility  of  5P4E  at 
compression  times  of  1  to  100  ps.  Although  not  clear  fi-om  the  latter  paper,  the  significance  of  a 
short  compression  time  is  that  the  process  approaches  the  adiabatic  idealization.  Therefore,  the 
short-time  compressibility  must  be  less  than  that  obtained  fi'om  long-time,  isothermal 
measurements  even  neglecting  structural  relaxation. 
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Table  3. 

Relative  Density  of  Water  by  Dilatometry 


P(P,T) 

p(o.T) 


0”C 

50"C 

95"C 

Pressure 

MPa 

DSD 

[25] 

DSD 

[26] 

[25] 

DSD 

[26] 

[25] 

0 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

1.000 

98 

1.046 

1.045 

1.040 

1.038 

1.039 

1.040 

1.041 

196 

1.080 

1.081 

1.072 

1.072 

1.076 

1.076 

294 

1.111 

1.112 

1.099 

1.100 

1.104 

1.105 

392 

1.136 

1.137 

1.121 

1.125 

1.126 

1.131 

490 

1.156 

1.159 

1.142 

1.143 

1.147 

1.153 

1.152 

1.154 

Table  4. 

Relative  Density  of  Santotrac  40 

P(P,T) 

p(o,T) 


Pressure 

MPa 

-15'’C 

25“C 

65“C 

0 

1.000 

1.000 

1.000 

100 

1.035 

1.043 

1.042 

200 

1.060 

1.074 

1.075 

300 

1.079 

1.097 

1.100 

400 

1.095 

1.113 

1.120 

500 

1.108 

1.132 

1.137 

600 

1.121 

1.151 

1.157 
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5.  NON-NEWTONIAN  RESPONSE 


The  search  for  the  relevant  constitutive  equations  which  relate  the  lubricant  stress  to  the 
flow  kinematics  in  concentrated  contact  has  occupied  the  interest  of  tribologists  working  in  EHD 
for  at  least  thirty  years.  The  goal  has  been  to  construct  experiments  to  verify  the  rheological 
models  and  to  provide  the  necessary  property  relations  so  that  the  complete  elastohydrodynamic 
solution  may  be  obtained.  Various  non-Newtonian  models  (see  [3]  for  a  review)  have  been 
advanced  which  provide  accurate  solutions  for  the  traction  over  some  operating  range  when  the 
required  rheological  properties  are  obtained  from  the  same  traction  data.  In  1972,  Dyson  [29] 
warned  that  EHD  traction  research  was  "enclosed  within  a  tight  circle"  of  fitting  parameters  to 
observations  without  consideration  for  measurements  made  outside  of  EHD.  Hopefully,  we  have 
broken  from  the  circle. 

We  may  now  generate  rheological  flow  curves  under  conditions  of  pressure,  temperature, 
and  rate  of  shear  which,  although  still  rather  restricted,  are  sufficiently  within  the  realm  of  EHD  to 
make  accurate  traction  predictions.  We  have  previously  interpreted  these  flow  curves  as  lubricant 
constitutive  behavior.  However,  in  light  of  the  recent  observation  [4]  of  mechanical  shear  bands 
operating  within  the  lubricant  film  concurrent  with  non-linear  shear  response,  this  interpretation 
must  be  accompanied  by  a  caveat:  Rheological  flow  curves  which  are  generated  in  plane  shear 
yield  an  empirical  rate  equation  which  is  useful  in  modeling  Couette  dominated  lubrication 
problems.  A  rigorous  analysis  of  the  EHD  problem  would  require  a  constitutive  equation  and  a 
slip  criterion  such  as  Mohr-Coulomb  [23]. 
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5.1  Viscous  Heating  in  Couette  Viscometry 

Presently,  the  most  useful  rheometer  configuration  for  investigation  of  high-pressure, 
high-stress  response  of  the  liquid  state  is  that  of  rotating  concentric  cylinders.  When  a  liquid  is 
sheared,  the  viscous  work  done  raises  the  temperature  of  the  liquid.  While  the  study  of  this 
phenomenon  is  of  itself  interesting  to  lubrication,  it  is  to  be  avoided  in  a  rheological  measurement 
because  constitutive  behavior  excludes  processes  which  result  fi-om  temperature  variation.  In 
previous  work  [30]  concerning  Couette  devices  the  authors  have  emphasized  the  importance  of  a 
low  Brinkman  No.  (through  primarily  a  small  shearing  gap)  to  mitigate  the  effect  of  the 
temperature  difference  within  the  liquid  film  and  fast  instrument  response  to  control  the 
temperature  of  the  surfaces  of  the  solid  boundaries  to  the  film.  An  alternate  experimental 
approach  is  to  perform  a  measurement  in  so  short  a  time  that  the  temperature  profile  in  the  film 
has  not  had  time  to  develop  but  at  a  late  enough  time  that  the  velocity  profile  is  fully  formed. 
Winter  [31]  showed  that  this  latter  approach  requires  that  the  Prandtl  No.  be  large. 

An  earlier  analysis  of  the  combined  effect  of  instrument  response  time  and  cylinder 
heating  [30]  modeled  the  cylinders  as  a  lumped  heat  capacity.  It  was  found  that  to  minimize 
errors  due  to  this  combination  the  measured  stress  history  should  show  no  thermal  softening  over 
a  period  of  time  equal  to  the  instrument  response  time.  Later,  the  evolution  of  the  radial 
temperature  profile  Avithin  the  cylinders  was  considered  [32].  A  detailed  example  of  a  numerical 
thermal  simulation  of  one  of  the  high-pressure  Couette  viscometers  used  in  this  laboratory 
follows. 

Only  conduction  in  the  radial  direction  is  considered.  An  appropriate  form  of  the  energy 
equation  for  the  metal  cylinders  is 
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(8) 


a^T  1  ar  _  p,c„  dr 

dr^  T  dr  k_  at 

n 

where  r  is  the  radial  coordinate,  t  is  time  and  p,  C,  and  k  are  the  density,  specific  heat  capacity  and 
thermal  conductivity,  respectively.  The  subscript,  m,  refers  to  the  metal  fi'om  which  the  cylinders 
are  fabricated.  Within  a  pressure  vessel  the  outer  surface  of  the  outer  cylinder  (or  cup)  is  in 
contact  with  low  conductivity  liquid.  Therefore  at  the  outermost  surface  the  boundary 
condition  is  set  at  dT/dr  =  0.  Symmetry  dictates  that  at  r  =  0,  dT/dr  =  0. 

The  shearing  gap  is  very  thin  compared  with  the  working  radius  so  that  curvature  of  the 
liquid  film  can  be  ignored.  A  preliminary  numerical  solution  for  the  liquid  film  with  isothermal 
boundaries  and  dissipation  showed  that  for  a  1  pm  film,  steady  state  temperature  distribution  was 
achieved  after  10  ps  of  shearing.  Therefore,  the  storage  term  in  the  energy  equation  for  the  liquid 
film  is  dropped.  Including  the  dissipation  term,  ty ,  we  have  for  the  liquid 

k0  +  Ty  =  O  (9) 

The  geometry  considered  is  shown  in  Figure  7.  The  liquid  being  sheared  is  in  the  circular 
gap  between  the  outer  and  inner  cylinders.  Only  the  outer  cylinder  rotates  with  a  velocity  which 
yields  a  specified  average  of  the  rate  of  shear,  y.  The  temperature  and  also  the  radial  heat  flux 
are  made  continuous  at  the  two  metal/liquid  interfaces.  Equations  (8)  and  (9)  in  their  respective 
regions  are  solved  numerically  with  the  above  boundary  conditions  and  with  initial  uniform 
temperature  (arbitrarily  set  to  zero  for  Figure  7)  to  obtain  temperature  distribution  and  shear 
stress  versus  time. 

The  result  of  a  sample  computation  is  shown  in  Figure  7  for  two  velocity  histories.  The 
velocity  of  the  outer  cylinder  was  increased  linearly  with  time  until  an  apparent  or  average  rate  of 
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shear  of  lO^s'l  was  reached  at  to  equal  to  10  and  50  ms.  The  liquid  rheology  was  Newtonian 
with  |x  =  10^  Pa  s  and  P  =  0.1°C‘^.  The  gap  was  1.0  pm.  The  liquid  and  metal  thermal 
conductivities  were  k  =  0.15  and  kjjj  =  225  W/m°C  and  for  the  metal  0^  =  3  x  10^  J/m^°C. 
The  temperature  profiles  within  the  metal  cylinders  are  shown  in  Figure  7  for  t  =  to-  The  liquid 
film  temperature  reached  a  maximum  very  near  the  surface  of  the  inner  cylinder.  The  reduction  in 
shear  stress  fi'om  isothermal  was  less  than  2%  for  to  =  10  ms  and  less  than  4%  for  to  =  50  ms. 
The  rotational  velocity  in  this  example  is  only  about  1  revolution  per  second  and  can  easily  be 
achieved  in  10  ms  by  a  stepper  motor.  Very  large  driving  torque  may  require  a  50  ms  linear 
velocity  ramp. 

In  early  high-pressure  Couette  devices  the  critical  response  time  was  often  that  of  the 
torque  transducer  which  provides  the  stress  measurement.  This  transducer  is  immersed  in  a  high- 
pressure  liquid  and  for  high  viscosity  the  transducer  responds  as  a  first  order  instrument  with  a 
time  constant,  tg,  which  is  proportional  to  damping  liquid  viscosity.  Very  fast  measurements  may 
be  obtained  by  accelerating  both  cylinders  and  then  arresting  the  rotation  of  one  cylinder  with  the 
torque  transducer  which  can  be  filled  with  a  low  ^^scosity  liquid.  Now,  the  velocity  response  is 
intimately  connected  with  the  transducer  response.  The  undamped  natural  frequency  for  the 
transducer  is  2.5  kHz.  Theoretically,  a  5  percent  settling  time  [33]  of  0.2  ms  can  be  achieved  for 
transducer  response  by  optimizing  the  damping  liquid  viscosity.  This  is  difficult  to  achieve  in 
practice  as  the  viscosity  of  the  damping  liquid  changes  with  the  test  temperature  and  pressure.  In 
any  event,  we  have  not  observed  an  improvement  in  measurements  typical  of  those  reported  here 
when  done  faster  than  50  ms.  The  temperature  increases  shown  in  Figure  7  present  no  problem  in 
interpretation  of  results. 
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Also,  by  not  considering  axial  conduction  in  the  cylinders  the  analysis  overestimates  the 
metal  temperature  rise.  Comparing  measured  stress  histories  for  a  Newtonian  liquid  with  the 
computed  stress  history  at  t  =  Is  we  found  the  actual  loss  of  shear  stress  compared  to  isothermal 
shear  was  only  64%  of  the  predicted  loss. 

The  temperature  distribution  of  Figure  7  was  incorporated  into  a  thermal-elastic  numerical 
analysis  of  each  cylinder  using  both  plane-strain  and  plane-stress  idealizations.  Thermal  expansion 
results  in  no  subtantial  change  in  the  operating  film  thickness  at  the  time  of  the  measurement. 

5.2  Some  Results  for  Pressures  to  300  MPa 

A  high-pressure,  high  shear  stress  Couette  viscometer  was  developed  [32]  in  this 
laboratory  to  operate  at  pressures  to  300  MPa.  The  steady  shear  response  of  various  liquid 
lubricants  [32]  and  a  grease  [34]  have  been  investigated.  For  the  grease,  an  empirical  stress 
equation  was  found  of  the  form 

=  xj  +  (gM-y)’  ,T  >  Xj,  (10) 

where  g  is  a  dimensionless  viscosity  enhancement  factor  approximately  equal  to  2,  q  is  a 
dimensionless  exponent  and  p  is  the  base  oil  viscosity  given  by  equation  (6).  For  q  =  0.5, 
equation  (10)  becomes  Casson's  equation  for  oil/particulate  mixtures. 

Flow  curves  are  presented  for  a  high-traction  fluid,  Santotrac  50  and  a  mineral  oil,  LVI 
260  in  Figures  8  and  9,  respectively.  Although  the  pressures  and  temperatures  are  comparable  the 
shear  response  is  quite  different.  The  traction  fluid  remains  Newtonian  to  a  shear  stress  of  at  least 
10  MPa  whereas  the  mineral  oil  responds  non-linearly  above  3  MPa.  Note  also  that  for  the 
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Figure  8.  Flow  Chart  for  Santotrac  50,  p  =  297  MPa. 


Figure  9.  Flow  Chart  for  Mineral  Oil,  LVI 260,  p  = 


mineral  oil,  the  stress  becomes  essentially  independent  of  rate  at  about  x  =  10  MPa  and  that  the 
temperature  dependence  diminishes  as  the  rate  dependence  diminishes  as  predicted  by  eqn.  (3). 
Because  of  their  extended  Newtonian  response  we  have  begun  to  use  the  cycloaliphatic  traction 
fluids  as  viscosity  standards  for  shear  stress  greater  than  1  MPa  in  calibrating  cylinder  gaps. 

5.3  A  New  High-Stress  Couette  Viscometer  for  600  MPa  Pressure 

An  increase  in  the  pressure  capability  of  previous  rotating  concentric  cylinder  viscometers 
is  necessary  to  investigate  the  shear  response  of  typical  lubricants  well  above  the  glass  transition 
temperature.  Translating  concentric  cylinder  rheometers  are  capable  of  pressures  above  1  GPa; 
however  the  rate  of  shear  is  so  low  that  the  viscosity  must  be  made  very  large  (by  operation  near 
Tg)  in  order  to  achieve  an  interesting  magnitude  of  shear  stress. 

We  recently  reported  [32]  the  development  of  a  Couette  viscometer  for  pressure  to  300  MPa 
and  the  results  of  the  previous  section  were  obtained  vdth  that  device.  The  pressure  capability  of 
the  Couette  technique  is  doubled  with  the  device  shown  in  section  view  in  Figure  10.  The  outer 
cylinder  (cup)  is  driven  by  an  external  stepper  motor  by  means  of  a  drive  shaft.  The  thrust  bearing 
which  prevents  expulsion  of  the  shaft  fi'om  the  vessel  has  been  moved  to  outside  of  the  vessel  so 
that  the  bearing  does  not  run  in  pressurized  liquid.  The  high-pressure  seal  is  now  a  simple  spring- 
loaded  packing.  There  is  no  need  for  high-pressure  electrical  feed-throughs  since  the  torque 
measurement  is  transmitted  out  of  the  vessel  optically.  The  inner  cylinder  (bob)  is  restrained  from 
rotation  by  a  torsion  bar  with  a  mirrored  surface  vrithin  a  glass  tube.  See  Figure  10. 

The  principle  of  the  torque  measurement  is  depicted  in  Figure  1 1 .  The  twist  of  the  torsion 
bar  results  in  a  deflection  of  a  laser  beam  which  is  detected  by  an  optical  position  sensor.  A  glass 
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Figure  11. 


Torque  Transducer  Operating  Principle  for  Rheometer  Shown  in 
Figure  10. 
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tube  (Figure  10)  provides  a  circular  interface  between  the  liquid  which  surrounds  the  mirror  and 
the  pressurizing  fluid  so  that  changes  in  refractive  index  of  individual  samples  need  not  be 
accounted  for.  For  highly  viscous  or  opaque  samples,  the  transducer  is  filled  with  a  low  viscosity 
fluorinated  oil  which  is  immiscible  with  the  sample. 


5.4  Results  for  Pressures  to  600  MPa 

A  flow  chart  is  presented  for  the  polyalphaolefin,  SHF  1001,  in  Figure  12.  The  response  is 
apparently  Newtonian  to  about  2  MPa  shear  stress,  although  the  effective  viscosity  below  2  MPa 
is  slightly  reduced  compared  to  falling  body  measurements. 

The  mineral  oil,  HVI  650,  has  been  the  subject  of  numerous  traction  experiments.  The 
shear  response  of  HVI  650  is  depicted  in  Figures  13  a  and  b.  The  curves  drawn  through  the  data 
points  for  23°C  and  44°C  represent  the  empirical  Carreau-Yasuda  equation 


x  =  iiY 


(11) 


Here,  a  is  a  dimensionless  parameter  which  controls  the  breadth  of  the  transition  from  Newtonian 
to  a  non-Newtonian  flow  regime  with  rate  sensitivity  coefficient  of  m.  For  m  =  0,  is  a  limiting 

shear  stress.  In  Figures  13  a  and  b,  we  used  a  =  1.5  at  23°C  and  a  =  2  at  44°C.  For  the  limiting 
stress,  the  form  proposed  by  Bezot,  et  al.  [35]  was  adopted 

'^L  =  Co  +  Cip  +  Cjp' 


(12) 
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Figure  12.  Flow  Chart  for  a  Polyalphaolefin. 


Figure  13a.  Flow  Chart  for  HVI 650. 
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Figure  13b.  Flow  Chart  for  HVI 650. 


Using  the  results  of  Figure  13b  and  Ref  [36]  and  [37],  we  obtained  cj  =  0.034  and  C2  =  2.1  x  10" 
^  MPa  with  Cq  arbitrarily  set  to  zero.  Note  that  for  HVI 650  at  44°C  (Figure  13a)  the  curves  are 
also  a  good  approximation  of  the  Sinh  Law, 

T  =  T„  sinh"‘([iY/x„)  (13) 

Where  is  a  stress  which  represents  the  limit  ofNewtonian  response  and  is  3.3,  5.2  and  6.5  MPa 
for  pressures  of 448,  517  and  586  MPa  respectively.  That  is,  /4  at  44°C. 

The  shear  response  of  the  traction  fluid,  Santotrac  40,  is  shown  in  Figure  14.  For  the 
pressure  of  414  MPa,  the  solid  curve  is  equation  (11)  and  the  two  broken  curves  are  equation 
(13)  for  of  6  and  25  MPa  respectively.  When  elastohydrodynamic  traction  curves  are 


interpreted  in  terms  of  equation  (13)  together  with  exponential  pressure- viscosity  the  traction 


must  be  suspect.  However,  if  the  pressure  of  a  rheological  measurement  is  chosen  carefully  it 
should  be  possible  to  obtain  agreement  between  the  Newtonian  limit  for  that  particular  pressure 
and  the  traction  gradient  since  the  Newtonian  limit  is  pressure  dependent  (eg.,  HVI  650  at  448 
MPa  and  44°C). 
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Figure  14.  Flow  Chart  for  Santotrac  40 


We  have  investigated  the  grease  from  Ref  [34]  and  its  base  oil  in  the  current  rheometer  to 
explore  departures  from  Newtonian  behavior  at  high  shear  stress.  In  Figure  15,  for  room 
temperature,  the  mineral  base  oil,  600P,  is  non-Newtonian  above  x  =  2  MPa  at  310  MPa  pressure. 
When  the  pressure  is  increased  to  517  MPa  this  oil  is  Newtonian  to  x  =  3  MPa.  The  CA7000, 
which  is  a  soap-thickened  grease  of  600P,  is  also  shown  in  Figure  15.  It  would  appear  from  these 
limited  data  that  the  shear  stress  for  the  grease  (at  high  stress)  may  be  obtained  by  multiplying  the 
base  oil  result  by  the  previously  [34]  obtained  viscosity  enhancement  factor,  g,  which  is  about  2. 
Clearly,  equation  (10)  is  not  applicable  here  since  it  reverts  to  Newtonian  at  high  shear  stress. 

The  High-Pressure  Impact  Viscometer  of  Wong,  et.  al.  [39]  entraps  a  quantity  of  liquid 
between  a  ball  and  a  plate.  Interferometry  is  used  to  determine  the  local  flow  rate  of  liquid 
leaking  from  the  entrapment  and  local  surface  distortion  which  yield  local  pressure  from  elasticity 
theory.  A  Rabinowitsch  Correction  for  slit  flow  gives  the  true  shear  rate  at  the  surface.  The 
authors  concluded  that  the  effective  viscosity  was  found  in  earlier  works  to  be  a  function  of  time 
because  the  shear  rate  varied  with  time  while  the  shear  stress  remained  at  the  limiting  value. 
Limiting  stress  type  behavior  was  observed  over  4  decades  of  shear  rate  with  m  »  0.01  for  LVI 
260. 

6.  EHD  TRACTION  CALCULATION 

We  should  expect  our  property  measurements  to  yield  reasonable  predictions  of 
concentrated  contact  traction.  Evans  [40]  generated  isothermal,  line-contact  traction  curves  for 
three  of  the  most  widely  investigated  liquid  lubricants:  5P4E,  Santotrac  50  and  HVI650.  These 
data  were  described  as  isothermal  since  the  disc  temperature  was  adjusted  to  provide  a  constant 
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Figure  15.  Flow  Chart  for  Grease  and  Base  Oil. 


estimated  average  film  temperature.  These  traction  data  are  presented  for  HVI650  in  Figure  16 
for  inlet  temperatures  of  40  and  60°C  and  average  pressures,  p,  of  0.47  and  0.63  GPa.  The 
sliding  velocity  is  AU  and  inlet  temperature  rise  is  3°C. 

For  the  traction  calculations  which  are  the  curves  in  Figure  16,  we  assumed  the  Hertzian 
pressure  distribution  and  integrated  the  Carreau-Yasuda  equation  (1 1)  across  the  contact  area  to 
obtain  the  average  shear  stress,  t.  The  dimensionless  parameter,  a,  was  specified  by  a  = 
63/(T+19°C)  which  yields  the  values  obtained  in  the  previous  section  without  approaching  the 
meaningless  condition  of  a  =  0  at  ordinary  temperatures.  Viscosity  and  limiting  stress  were 
obtained  from  equations  (6)  and  (12).  Results  were  insensitive  to  the  selection  of  the  rate 
sensitivity,  m,  from  0  to  0.03.  Although  the  predicted  curve  at  40°C  rolls  over  more  quickly  than 
the  measured  traction,  the  general  agreement  is  good.  Considering  the  great  differences  between 
the  two  techniques  and  the  assumptions  involved  (eg.,  Hertzian  pressure)  it  might  be  unreasonable 
to  expect  to  do  better. 

7.  NEW  OBSERVATIONS  OF  MECHANICAL  SHEAR  BANDS  IN  SIMPLE  PLANE 
SHEAR 

The  High-Pressure  Flow  Visualization  Cell  which  we  use  to  optically  observe  shear  bands 
has  been  described  in  detail  in  Ref  [4]  as  well  as  the  experimental  procedures  and  preliminary 
results.  It  was  shown  that  at  a  critical  shear  stress  which  marked  the  limit  of  linear  response,  an 
intermittent  slip  occurred  across  visible  bands.  The  existence  of  two  types  of  bands  with  differing 
orientation  is  now  recognized  [41].  Figure  17  illustrates  the  orientations  of  the  two  types.  New 
observations  which  give  insight  into  the  slip  mechanism  are  presented  in  the  next  sections. 
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Figure  16.  Measured  and  Predicted  Traction  Curves  for  HVI 650. 
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Figure  17.  Two  types  of  mechanically  induced  shear  bands. 


7.1  Birefringence 

Flow  birefringence  has  been  useful  for  stress  analysis  in  flowing  liquids  in  a  manner 
analogous  with  photoelasticity  in  transparent  solids  [42].  Molecular  orientation  under  stress  gives 
rise  to  optical  anisotropy  -  the  refractive  index,  n,  is  different  in  different  directions.  Birefringence 
observations  can  be  used  to  determine  the  principal  normal  stress  difference,  Aa,  and  the 
orientation  of  the  principal  stress  axes  if  a  suitable  calibration  of  the  birefringence,  An,  can  be 
obtained.  Such  a  calibration  is  not  possible  with  the  flow  cell.  However,  some  generalizations 
regarding  the  stress  distribution  in  the  flow  field  may  be  obtained. 

The  High-Pressure  Flow  Visualization  Cell  was  placed  between  crossed  polarizers.  White 
light  was  filtered  to  produce  a  narrow  band  at  605  nm.  The  polyphenyl  ether,  5P4E,  was  the 
model  lubricant.  Temperature  was  23°C  and  pressure  was  241  MPa.  The  shearing  force  was 
increased  with  time  until  shear  bands  were  observed.  The  interference  pattern  shown  in  Figure  18 
was  obtained  just  before  the  shear  localized.  This  fringe  pattern  disappeared  at  the  time  the  shear 
bands  formed.  Similar  experiments  were  performed  with  the  crossed  polarizers  placed  at  eight 
different  angles  with  respect  to  the  shear  direction  equally  spaced  over  a  span  of  180°  in  an 
attempt  to  determine  the  extinction  angle.  No  effect  of  orientation  of  the  polarizer/analyzer  pair 
was  noted  on  the  contrast  or  geometry  of  the  fringes  or  shear  bands. 

Apparently,  the  high-pressure  windows  which  are  sapphire  are  performing  as  fractional 
wave  plates  to  eliminate  the  extinction  angle  effect  [42]  and  so  it  was  not  possible  to  obtain  the 
principal  stress  directions.  It  is  possible  to  obtain  the  stress  distribution  in  a  qualitative  sense.  If 
the  stress-optical  function  is  monotonic  [43],  then  the  fiinge  order  is  a  representation  of  the 
relative  magnitude  of  the  principal  normal  stress  difference  (and  the  principal  shear  stress). 
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Referring  again  to  Figure  18,  the  zeroeth  order  fringe  occurs  in  the  liquid  reservoir  far  to  the  right 
in  the  figure.  The  principal  shear  stress  should  be  a  maximum  within  the  highest  order  fiinge  in 
the  region  marked  A.  Notice  that  the  shear  stress  is  uniform  in  the  shearing  gap  at  positions 
farther  than  about  3  times  the  film  thickness  into  the  gap. 

It  was  previously  noted  [4]  that  the  first  band  nucleates  at  the  point  A  and  runs  to  the 
point  marked  B  in  Figure  18.  Two  conclusions  may  be  drawn  from  this  observation.  1)  The 
reason  for  the  observation  of  the  first  bands  near  the  entrance  is  that  the  shear  stress  is  greatest 
there.  2)  A  band  which  nucleates  within  a  region,  A  of  locally  high  stress  (compared  to  the 
average  stress  within  the  field)  will  continue  to  run  through  a  region  of  locally  lower  (than  the 
average)  stress.  The  significance  of  the  latter  is  important.  Once  initiated  the  shear  band  defines  a 
plane  on  which  slip  is  more  easily  accommodated  than  in  the  bulk  of  the  material  -  it  becomes  a 
“weak  spot”  and  the  slip  is  arrested  only  by  pinning  at  the  solid  boundaries  of  the  film. 

7.2  Persistence  of  Shear  Bands 

When  mechanically  induced  shear  bands  are  studied  in  solid  amorphous  polymers,  the 
observations  are  usually  made  on  sectioned  specimens  [24]  long  after  the  deformation  responsible 
for  the  bands  has  ceased.  This  is  possible  because  the  image  of  the  band,  which  is  probably  a 
result  of  damage  or  dilatation,  persists.  In  the  liquid  lubricants  this  persistence  time  is  shorter.  In 
Figure  19  we  have  plotted  the  time  for  which  a  shear  band  was  visible  in  5P4E  after  shearing 
(persistence  time)  versus  the  viscosity  at  test  conditions.  If  we  define  the  mechanical  shear 
relaxation  time  as  the  ratio  of  viscosity,  p,  to  shear  modulus,  G,  then  the  straight  line  plotted  in 
Figure  19  is  the  mechanical  relaxation  time  for  G  =  1  GPa.  We  may  define  a  characteristic  time 
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Figure  19.  Time  for  which  a  shear  band  is  visible. 
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for  thermal  difiiision  as  h^/D  where  h  is  the  thickness  of  the  band  and  D  is  thermal  dififiisivity.  D 
is  essentially  independent  of  p  and  is  typically  about  10'^  m^/s  for  liquid  lubricants.  If  a  shear  band 
is  at  most  2pm  thick,  then  the  characteristic  time  for  thermal  difiiision  is  at  most  40ps.  The 
correlation  of  persistence  (of  the  order  of  seconds)  vwth  mechanical  relaxation  is  much  better  than 
for  thermal  difiiision.  Clearly  the  image  of  the  shear  band  is  a  result  of  mechanical  damage  or 
dilatation  and  not  temperature. 


is  1.63.  Inserting  these  values  into  eqn.  (16)  we  find  that  the  rate  of  change  of  n  with  respect  to  e 
is  -0.3  for  143AD  and  -0.8  for  5P4E.  So  for  the  same  dilatation  the  refractive  index  changes 
three  times  less  for  the  perfluoropolyalkyl  ether  than  it  does  for  the  polyphenyl  ether.  Thus  the 
relative  degree  of  contrast  in  shear  band  images  is  consistent  with  an  argument  that  a  persistent 
dilatation  has  occurred  along  the  band. 
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7.4  Surface  Roughness  Effect 

It  might  be  argued  that  the  roughness  of  the  solid  boundary  provides  a  nucleation  site  for 
shear  bands.  The  moving  shaft  which  performs  as  one  (the  moving  one)  of  the  solid  boundaries 
has  two  sides  -  each  of  which  may  be  utilized  in  an  experiment.  One  side  was  left  as  machined 
with  an  rms  roughness  of  1.0  pm.  The  other  side  was  polished  to  a  roughness  of  0.03  pm.  The 
stationary  surface  has  a  roughness  of  0.3  pm  rms.  Both  sides  of  the  moving  shaft  has  been 
utilized.  No  difference  in  the  character  of  the  bands  or  the  manner  in  which  they  developed  was 
noted.  We,  therefore,  cannot  associate  surface  roughness,  at  least  up  to  1  pm,  with  shear 
localization  in  liquid  films. 

8.  SLIP  CRITERION 

8,1  Analogy  to  Glassy  Polymers 

In  previous  work  [4],  we  have  associated  mechanically  induced  shear  bands  with  the  rate- 
independent  shear  stress  which  is  observed  in  liquid  lubricants  under  pressure  in  disc  machines  and 
rheometers.  As  a  critical  stress  which  is  roughly  proportional  to  pressure  is  approached,  an 
increasing  amount  of  the  deformation  in  an  otherwise  rate-dependent  matrix  is  accommodated  by 
intermittent  slip  along  inclined  shear  planes.  When  this  critical  (limiting)  stress  is  reached,  any 
increase  of  the  apparent  rate  of  shear  can  be  accommodated  by  an  increased  production  of  shear 
bands  without  changing  the  stress:  the  response  is  rate-independent. 

In  practice,  for  glassy  polymers,  the  strain  rate  sensitivity  coefficient  m  =  ,  and 

be  quite  low  when  shear  bands  are  operating.  Here,  x  is  shear  stress  and  y  is  rate  of  shear. 
G’Sell  [44]  reports  that  m  =  0.027  for  polycarbonate  with  shear  bands.  Such  a  low  rate  sensitivity 
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(6%  per  decade)  would  possibly  be  undetectable  in  high-pressure  rheometers  and  disc  machines; 
hence  the  limiting  shear  stress  idealization. 

Glassy  polymers  are  known  to  exhibit  a  pressure  dependent  shear  strength  which  at  high 
pressure  follows  the  rule  [24] 

Xy  =  Xyo  +  ^yP  (17) 

where  Xy  varies  from  0.1  to  0.25  and  p  is  pressure.  A  similar  rule  is  known  to  apply  to  the 
limiting  shear  stress  of  liquid  lubricants,  although  with  a  lower  proportionality  constant. 

8.2  Mohr-Coulomb  Criterion 

The  Mohr-Coulomb  criterion  predicts  that  slip  will  occur  along  any  plane  on  which  the 
ratio  of  shear  stress,  Xe,  to  compressive  normal  stress,  -oe,  attains  the  magnitude  of  a  material 
friction  coefficient,  t].  Referring  to  Figures  20  and  21,  xe  and  oe  may  be  resolved  from  the  shear 
stress,  T,  and  normal  stresses  Ox  and  Oy  which  are  oriented  along  and  perpendicular  to  the  solid 
boundaries  in  a  plane  shear  experiment.  The  mean  mechanical  pressure,  p,  is  defined  by  p  =  -!4(ax 
+  Oy)  and  X  =  x/p. 

In  general,  the  normal  stresses  will  not  necessarily  be  equal  and  we  will  quantify  the  first 
normal  stress  difference,  Ni  =  Ox  -  Oy  by  ^  =  Ni/2p.  The  parameters  X  and  ^  are  dimensionless. 
We  may  resolve 


-Oe  =  p  +  ^pcos20  -  xsin20 

(18) 

Xe  =  xcos20  +  ^psin20 

(19) 

The  slip  criterion  can  be  represented  graphically  on  the  Mohr’s  circle  plot  in  Figure  20. 
The  Mohr’s  circle  of  stress  must  fit  within  the  envelope  defined  by  the  broken  lines.  These  lines 
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Figure  20.  Mohr’s  circle  representation  of  slip  criterion. 


Figure  21.  Definitions  of  angles  and  stresses. 


must  curve  away  from  the  a  axis  near  p  =  0  to  allow  a  non-zero  shear  stress  there.  The  included 
angle  of  the  envelope  is  twice  the  material  friction  coefiBcient  "n.  Now,  slip  will  occur  when  the 
circle  of  stress  increases  to  tangency  with  the  envelope  and  the  orientation  of  the  band  will  be 
such  that  |Xe  /  Gel  is  maximized.  This  is  satisfied  by  the  two  radii  of  the  circle  drawn  to  the  points 


of  tangency.  These  radii  represent  the  two  shear  band  angles,  0i  and  02  in  Figure  21. 
Analytically,  the  shear  band  angles  can  be  found  by  setting 


d  1 


d0 1  -c, 


=  0 


(20) 


8=e. 


which  leads  to 


X?  +  -  Xsin20^  +  ^cos20g  =  0 


(21) 


The  two  solutions  for  0c  are  the  shear  band  angles  0i  and  02  for  the  fifst  and  second  types 
respectively.  The  material  fiiction  coefficient  is 


X  cos  20,  +  Esin  20, 
n  = - 2 — - £— 

1  +  E cos 20,  -  Xsin20, 


(22) 


If  the  first  normal  stress  difference  Ni  were  zero,  then  the  two  shear  band  angles  would  be 
complimentary  angles  (sum  to  90°)  defined  by  solutions  of  0,  =  j^sin"*  X .  Clearly  the  shear 


band  angles  in  Figure  17  are  not  complimentary.  It  will  be  necessary  to  invoke  the  first  normal 
stress  difference  to  reconcile  theory  with  observation. 

For  example,  from  Figure  17  for  the  polyphenyl  ether,  0i  =  19°  and  02  =  103°.  Then  X  = 
0.089  and  E  =  0.056  are  obtained  from  equation  (21).  This  value  of  X  is  in  agreement  with  the 
measured  [32]  limiting  shear  stress  if  X  is  applied  to  a  linear  equation  like  (17).  The  first  normal 
stress  difference  obtained  is  a  little  more  than  the  shear  stress.  This  may  seem  excessive,  since 


55 


consideration  of  the  normal  stress  difference  has  not  been  necessary  to  predict  bearing  load 
capacity;  however.  Tanner  [45]  showed  that  for  typical  lubrication  flows  the  fractional  increase  in 
load  capacity  is  unaffected  by  Ni  of  the  order  of  x.  The  experimental  measurement  of  Ni  at  high 
pressure  is  challenging  and  has  not  previously  been  attempted  -  possibly  because  of  the  expected 
irrelevance  to  bearing  load  capacity.  Measurement  of  Ni  could  confirm  the  applicability  of  Mohr- 
Coulomb  and  might  be  considered  a  research  priority. 

In  the  example  above,  the  material  fiiction  coeflBcient,  r\  =  0.106.  The  value  of  t)  is  simply 
the  difference  between  nil  and  the  included  angle  between  the  two  types  of  band  (in  radians). 
Since  this  included  angle  is  very  nearly  7t/2,  the  determination  of  (and  X)  from  observations  of 
bands  is  sensitive  to  the  accuracy  of  the  shear  band  angle  measurement. 

The  first  normal  stress  difference,  Ni,  in  simple  shear  was  the  subject  of  conjecture  and 
controversy  early  in  the  study  of  high-pressure  rheology  of  lubricants;  however  there  has  been  no 
attempt  to  measure  this  effect  under  pressure.  Gleissle  [46]  empirically  associated  flow 
instabilities  with  a  critical  ratio  of  1st  normal  stress  difference  to  shear  stress.  He  found  that  for 
melt  fracture  in  capillaries  Ni/x  =  4.6  and  alluded  to  an  unpublished  analysis  which  predicted 
instability  for  Ni/x  >  4.  Abdel-Khalik,  et  al.  [47]  present  a  method  of  predicting  Ni  fi-om 
measured  flow  curves.  They  considered  (among  others)  eqn  (11)  with  a  =  2  and  m  =  0.  Their 
result  is  that  Ni/x  =  1.26  for  [IY/Xl  =  0.35  which  is  somewhat  consistent  with  the  appearance  of 
bands  at  the  onset  of  non-linearity.  For  Ni/x  =  4.6,  Gleissle’s  limit  for  stable  flow  in  capillaries, 
tty  /  Xl  =  1.2 .  In  a  purely  analytical  study  [48]  of  a  constitutive  law  which  combines  the  Maxwell 
fluid  model  and  a  compressible  plasticity  model  we  have  shown  that  the  shear  band  orientation  as 
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determined  by  the  method  of  characteristics  is  independent  of  pressure  and  consistent  with  the 
values  observed  experimentally.  (See  Appendbc) 

9.  RHEOLOGICAL  MODELS 

In  previous  work  [3]  the  authors  showed  from  experimental  measurements  that  the 
Maxwell  Model,  which  sums  elastic  and  viscous  strain  components,  correctly  describes  the 
transient  liquid  response  under  pressure.  The  simplest  incompressible  form  of  the  Maxwell  model 
requires  two  rheological  properties  -  for  example:  the  limiting  elastic  shear  modulus,  G,  and  the 
limiting  low  shear  viscosity,  p. 

A  full  rheological  model  for  elastohydrodynamic  lubrication  was  developed  from 
observations  of  the  liquid  lubricant  response  observed  in  high-pressure  rheometers.  Here,  dij  is 
the  rate  of  deformation  tensor,  xij  is  the  deviatoric  stress  tensor  and  Te  is  the  von  Mises  stress.  In 
formulating  this  equation  we  adopted  the  Stokes’  Condition  -  the  mean  mechanical  pressure,  p, 
being  set  equal  to  the  thermodynamic  pressure,  pr.  However,  in  writing  the  full  model  in  Ref  [3] 
we  also  set  the  second  coefficient  of  viscosity  equal  to  zero  and  this  is  inconsistent  with  Stokes’ 
Condition.  We  suggest  tentatively  that  a  term  be  added  to  the  right-hand  side  of  the  equation  so 
that  the  full  model  now  reads 


+ 


2p 


F(T,p,xJ 


+ 


2p 


(Pt  -  p) 


(23) 
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where  K  is  now  the  bulk  modulus  of  the  glass  and  pr  is  that  pressure  which  yields  the 
instantaneous  density  of  the  liquid  when  used  in  the  equilibrium  equation  of  state.  Stoke’ s 
Condition  is  now  abandoned  and  compressional  viscoelasticity  is  addressed  explicitly.  The  second 
coefficient  of  viscosity  is  zero  which  sets  the  bulk  viscosity  equal  to  2/3  p.  This  choice  is  made  so 
that  the  short  time  isothermal  compressibility  is  1/K  and  the  long  time  compressibility  is  obtained 
fi'om  the  state  equation.  The  function  F(T,p,  te)  is  an  empirical  rate  relation. 

Equation  (23)  represents  unfinished  business  in  two  respects.  A  first  normal  stress 
difference  is  not  explicitly  introduced.  The  use  of  the  Jaumann  time  derivative  in  (23)  will  result 
in  Ni/x  of  the  order  of  x/G  if  F  =  1  [49].  Thus  the  magnitude  of  shear  stress,  x,  must  be  close  to 
that  of  the  shear  modulus,  G,  to  yield  the  first  normal  stress  difference  required  by  the  Mohr- 
Coulomb  theory  and  observed  shear  band  angles.  This  is  too  great  a  value  of  x  since  Xi/G  is 
approximately  1/30.  Secondly,  we  now  know  that  the  non-linear  (in  Xe)  form  of  F(T,p,  Xe)  is  at 
least  in  part  due  to  shear  localization.  Constitutive  behavior,  by  definition,  excludes  the  effect  of 
localization.  For  many  simple,  low  molecular  weight  base  stocks  it  may  be  most  correct  to  set  F  = 
1  and  apply  a  slip  criterion  such  as  Mohr-Coulomb.  Then  the  non-linear  behavior  which  was 
introduced  previously  through  F  is  a  consequence  of  the  distribution  of  t|  through  the  material  and 
the  slip  velocity. 

For  simple  cases  where  the  flow  is  steady  simple  shear,  an  empirical  stress  equation  such 
as  we  have  advanced 

T  =  T,(l-e-''''-)  (24) 

is  sufficient.  Here  the  rate  sensitivity  coefficient,  m,  goes  to  zero  as  py/  Tl  becomes  large.  For 
solid  polymers  where  shear  bands  are  operating,  m  is  small  but  not  zero.  Also,  recent 
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measurements  [32]  indicate  that  the  transition  from  Newtonian  to  “rate-independent”  behavior  for 
high  molecular  weight  lubricants  is  broader  in  shear  rate  than  that  which  is  described  by  equation 
(24).  These  deficiencies  may  be  removed  by  using  the  Carreau-Yasuda  form  introduced 
previously 


The  dimensionless  exponent,  a,  controls  the  breadth  of  the  transition  and  rate  sensitivity 
coefiBcient,  m,  appears  explicitly.  Note  that  equation  (11)  is  equivalent  to  the  Elsharkawy  and 
Hamrock  [50]  model  when  m  =  0.  For  m  very  small  (-0.01)  but  not  zero,  equation  (11)  fits 
experimental  results  well  while  removing  the  singularity  which  causes  problems  in  numerical 
simulations. 


10.  MECHANICAL  SHEAR  BANDS  IN  NON-VISCOMETRIC  FLOW 
10.1  Shear  Bands  in  Wedge  Flow 

Visualization  of  shear  bands  has  to  date  been  accomplished  between  parallel  plates.  The 
interpretation  of  the  flow  field  and  stress  state  leading  up  to  localization  is  the  least  ambiguous  for 
plane  Couette  shear.  This  is  not  however  the  general  case  for  lubrication  flows.  To  investigate 
the  generation  of  shear  bands  in  wedge  flow  we  fabricated  a  new  stationary  shaft  for  the  High- 
Pressure  Flow  Visualization  Cell.  The  stationary  surface  (the  lower  surface  of  Figure  22)  was 
ground  at  an  angle  of  5°  with  the  direction  of  motion  of  the  moving  (upper  of  Figure  22)  surface. 
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Figure  22. 


5P4E.  138  MPa,  S^C 


Mechanical  Shear  Bands  between  Converging  Plates  with  Sliding.  Lower 
Tilted  Surface  is  Stationary. 
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Observation  of  birefringence  during  shear  through  crossed  polarizers  clearly  showed  a  stress 
gradient  across  the  film  as  expected.  Shear  bands  appeared  as  shown  in  Figure  22  for  SP4E. 
These  bands  are  similar  to  those  observed  between  parallel  plates.  However,  the  analysis  of  the 
previous  section  can  not  be  applied  because  of  the  non-uniform  stress  field  in  the  film. 

10.2  Extra-High-Pressure  Flow  Visualization  Cell 

The  simple  plane  shear  configuration  is  preferable  for  the  interpretation  of  flow 
visualization.  However,  the  pressure  capability  is  limited  by  the  fiiction  of  shaft  seals.  A  new 
Extra-High-Pressure  Flow  Visualization  Cell  is  shown  in  section  view  in  Figure  23.  Two 
windows  are  provided  -  one  for  illumination  as  shown  and  one  for  viewing.  The  entire  liquid 
sample  is  compressed  by  the  motion  of  the  piston  and  rod  to  generate  the  ambient  pressure  of  up 
to  1  GPa.  A  hydraulic  actuator  pushes  against  the  rod.  Pressure  measurement  is  obtained  from 
the  force  on  the  rod  and  a  simple  hysteresis  calibration  of  seal  fiiction.  A  pin  (1.2  mm  diameter) 
is  attached  to  the  piston  so  that  after  some  compression  of  the  liquid,  a  squeeze  film  is  formed 
between  the  pin  and  the  anvil.  The  ambient  pressure  just  before  squeeze  can  only  be  adjusted  by 
varying  the  initial  liquid  volume.  The  illumination  window  is  tilted  so  that  due  to  refraction  light 
will  approach  the  gap  between  the  pin  and  anvil  at  a  glancing  angle. 

10.3  Shear  Band  Observations  in  Squeeze  Film 

When  a  mineral  oil  (N1  from  [18,  37])  is  pressurized  to  1.0  GPa  at  room  temperature  and 
then  squeezed  between  the  pin  and  anvil,  shear  bands  develop  near  the  edge  of  the  gap  as  shown 
at  the  top  of  Figure  24.  The  dashed  white  lines  complete  the  outlines  of  the  pin  and  the  anvil 
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gh-Pressure  Flow  Visualization  Cell. 
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below  the  pin.  The  axis  of  the  pin  is  vertical  and  to  the  left  in  the  upper  picture.  For  comparison, 
shear  bands  in  solid  polymethylmethacrylate  sectioned  from  a  specimen  deformed  in  plane-strain 
compression  are  shown  in  the  bottom  of  Figure  24  (taken  from  [24]).  The  similarity  in 
appearance  and  orientation  is  striking.  The  glass  transition  pressure  at  room  temperature  is  0.6 
GPa  for  N1  so  the  comments  of  Section  4  regarding  hydrostatic  pressure  apply.  The  reason  that 
the  bands  are  not  at  90°  to  one  another  (as  are  the  principal  shear  directions)  for  both  the  N1  and 
the  polymer  has  been  attributed  [24]  to  dilatation  or  the  pressure  dependence  of  yield  stress.  The 
slip  planes  are  turned  in  a  direction  of  reduced  compressive  normal  stress  for  which  slip  is  more 
favorable.  Shear  bands  were  also  observed  in  HVI  650  at  0.8  GPa  in  the  present  device.  These 
results  are  relevant  to  squeeze  film  lubrication  and  plate  impact  and  diamond  anvil  experiments. 

11.  ADIABATIC  SHEAR  BANDS 

In  previous  sections  we  have  discussed  mechanical  shear  bands.  A  very  different  type  of 
shear  localization  is  the  subject  of  this  section.  We  present  the  experimental  observation  and 
numerical  simulation  of  a  thermally  driven  localization  similar  to  the  “adiabatic  shear  bands”  found 
during  high  rate  deformation  of  metals  [51].  We  have  adopted  the  term  “adiabatic”  from  the 
usage  in  plasticity  literature,  although  the  term  is  not  strictly  correct. 

11.1  Background 

In  1967,  Plint  [52]  recognized  that  the  inhomogeneous  shear  in  an  elastohydrodynamic 
film  which  exists  due  to  viscous  heating  might  localize  at  or  near  the  mid-plane  where  the 
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MINERAL  OIL  In  1GPa  Ambient 


PMMA  POLYMER  in  1atm  Ambient 


Figure  24  Shear  Bands  in  Mineral  Oil  Squeezed  between 
Pin  and  Anvil  (Top)  and  Shear  Bands  in 
Solid  Polymer  in  Compression 
(Bottom  from  [24]). 
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temperature  is  greatest.  If  the  shear  stress,  t,  is  uniform  across  the  film,  the  viscous  dissipation 
will  be  approximately  exponential  with  temperature.  The  lubricating  film  is  generally  a  poor 
thermal  conductor.  Therefore,  as  the  temperature,  T,  rises  locally  the  dissipation  which  generates 
the  temperature  field  rises  rapidly  with  the  local  rate  of  shear,  y ,  and  conditions  favor  an  unstable 


temperature  profile. 

The  parameter  which  characterizes  the  thermal  behavior  of  liquid  films  in  plane  Couette 
shear  is  the  dimensionless  Brinkman  number  (also  known  as  GrifBth  or  Nahme  number) 


where  h  is  the  thickness  of  the  film,  po  is  the  initial  viscosity  at  T  =  To,  k  is  the  thermal 
conductivity  of  the  liquid  and  P  is  the  temperature-viscosity  coeflScient. 

It  has  been  long  recognized  in  the  rheology  literature  that  a  low  value  of  Br  is  a  necessary 
condition  to  avoid  the  influence  of  viscous  heating  in  continuous  shear  viscosity  measurement 
[53],  However,  the  significance  of  the  Brinkman  number  is  not  well  recognized  in  tribology.  The 
Brinkman  number  appears  as  the  coefficient  of  the  dissipation  term  in  the  dimensionless  energy 
equations  [12,  45].  A  steady  state  solution  to  the  energy  equation  exists  for  Br  <  3.51  when  the 
temperature-viscosity  function  is  exponential  and  for  Br  <  7t^  when  the  viscosity  varies  linearly 
with  temperature  [31].  Experimental  investigations  of  lubricant  constitutive  behavior  for  EHD 
conditions  which  may  be  regarded  as  isothermal  are  performed  at  Br  <  .3  in  high  pressure 
rheometers  [30,  32]  and  disc  machines  [36].  Thermal  localization  phenomenon  are  expected  to 
manifest  at  large  (>3.51)  values  ofBr. 
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11.2  Experiment 

The  High-Pressure  Flow  Visualization  Cell  has  been  described  in  detail  elsewhere  [4], 
This  device  is  shown  schematically  in  Figure  25.  The  liquid  lubricant  sample  is  sheared  in  the  gap 
between  a  flat  on  a  moving  shaft  and  the  end  of  a  stationary  shaft.  The  gap  is  illuminated  and 
viewed  in  a  direction  perpendicular  to  the  direction  of  motion  and  parallel  to  the  plane  of  the 
bounding  surfaces.  The  moving  shaft  is  pushed/pulled  by  a  gas  actuated  piston  and  this  shaft  must 
pass  through  high-pressure  closure  seals.  Motion  is  not  continuous  reciprocation,  but 
unidirectional  for  a  given  run. 

For  a  shear  rheological  measurement  during  which  rheological  properties  vary  with  time, 
one  would  choose  to  design  the  experiment  following  one  of  two  idealizations  [45]:  1.  rate  (of 
strain)  control  where  the  relative  velocity  is  kept  constant,  or  2.  force  (stress)  control  where  the 
average  shear  stress  is  imposed  as  constant.  Whereas,  it  is  known  that  rate  control  will  not  yield 
thermal  runaway  [45],  we  intend  to  show  that  in  a  stress-controlled  experiment  localization  can 
indeed  take  place.  The  device  shown  (a  constant  actuating  gas  pressure  applied  to  the  piston) 
would  be  ideal  for  force  controlled  experiments;  however,  there  exists  some  fiiction  in  the  high- 
pressure  seals,  which  is  known  to  have  a  viscous  component. 

The  liquid  lubricant  chosen  for  this  study  is  the  five-ring  polyphenyl  ether,  5P4E  for  which 
there  is  probably  the  most  high  pressure  property  data  available.  The  thermal  properties,  k=0.15 
W/m°C  and  pC=1.85  MJ/m^  ”C,  of  Richmond,  et  al.  [54]  and  the  rheological  properties  of 
Yasutomi,  et  al.  [18]  were  used.  Refractive  index  was  measured  in  a  high  pressure  refractometer 
described  by  Bair  [13].  In  addition,  the  properties  of  this  material  have  a  great  sensitivity  to 
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Figure  25.  Flow  Visualization  Configuration.  (Ulumination  Through  Lower 
Window.) 


pressure  so  that,  at  relatively  low  pressure,  behavior  is  observed  which  is  representative  of  many 
materials  at  much  greater  pressure. 

Regulated  compressed  nitrogen  gas  was  applied  to  the  left  side  of  the  piston  in  Figure  25 
with  sufficient  pressure  (3  MPa)  to  develop  an  initial  velocity  of  0.75  mm/s.  The  sample  pressure 
was  172  MPa  and  the  initial  temperature  22°C  yielding  a  viscosity  |Jo  =  3  MPa  s  and  3  =  0.6°C'\ 
Although  the  piston  force  was  constant,  the  force  contributed  by  the  seals  increased  with  sliding 
velocity  so  that  the  average  stress  applied  at  the  shearing  gap  decreased  with  velocity  (and  time). 
A  standard  (30  frames/second)  video  camera  recorded  the  video  prints  shown  in  time  sequence  in 
Figure  26.  The  lower  black  object  is  the  stationary  boundary  and  the  upper  is  the  moving  (left  to 
right)  boundary.  The  gap  or  clearance  is  150  pm  ,  and  x  is  obtained  fi'om  the  product  of  initial 
viscosity  and  the  initial  rate  of  shear,  to  give  Br  =  6.75.  The  central  one-third  of  the  length  (in  the 
flow  direction)  is  captured  in  each  print. 

The  velocity  history  is  shown  by  the  points  in  Figure  27  -  obtained  by  differentiating  the 
signal  from  an  LVDT  on  the  moving  shaft.  Past  210  ms  the  velocity  exceeded  the  measurement 
capability. 

A  striking  feature  in  the  sequence  of  images  in  Figure  26  is  the  abrupt  appearance  of  a 
dark  band  within  the  liquid  film  at  260  ms.  In  the  previous  fi'ame  at  230  ms  the  band  is  absent. 
However  the  interior  of  the  film  has  begun  to  darken.  In  less  than  30  ms  a  dark  band  has 
appeared.  This  feature  faded  and  disappeared  within  10s  after  the  interruption  of  shear.  We  show 
in  the  next  section  that  the  dark  band  is  consistent  with  the  appearance  of  a  hot  liquid  layer 
surrounded  by  cooler  liquid.  Light  transmitted  through  the  film  should  be  refi'acted  from  a  hot 
region  to  cool  regions.  Effectively,  two  cool,  viscous  liquid  films  are  sliding  against  a  hot  low 


68 


# 


Figure  26, 


5P4E,22C  INITIAL  TEMP,  172  MPa 

Visualization  of  an  Adiabatic  Shear  Band.  Time  is  Elapsed  from  Start 
of  Shear. 
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viscosity  film.  The  curvature  of  the  band  is  apparently  due  to  cool  liquid  being  drawn  into  the 
shear  gap  by  the  moving  shaft  at  the  upper  left  of  the  prints  in  Figure  26.  That  the  band  which 
suddenly  appears  represents  shear  localization  will  be  shown  by  the  analysis  which  follows. 


11.3  Thermal  Analysis  of  Plane  Couette  Shear  Under  Stress-Controlled  Conditions 

The  statement  of  the  conservation  of  energy  translates  into  the  following  commonly 
known  as  the  energy  equation 

P‘:,^={kT4  +  <|.,  (26) 

where  the  left  member  of  the  above  equation  represents  the  energy  transport  by  convection  with 
D/Dt  representing  the  material  derivative  and  pcp,  the  volume  heat  capacity.  The  first  term  on  the 
right  denotes  the  energy  transfer  by  conduction  followed  by  the  viscous  dissipation  term  which  is 
<|)  =  try. 


For  the  geometry  under  consideration,  the  energy  equation  reduces  to 


aT  ,  a^T  ^ 


(27) 


where  y  is  the  coordinate  normal  to  the  direction  of  motion  in  the  direction  of  the  film  thickness. 


If  the  Newtonian  assumption  is  adopted,  the  viscous  dissipation  term  becomes  (b  =  — , 


where  the  viscosity,  p,  varies  exponentially  with  temperature  (for  small  T  -  To)  prescribed  as 


-P(T-T,) 


(28) 


The  dimensionless  form  of  the  equation  (27)  for  a  stress-controlled  experiment  takes  on  the  form 
given  below 
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(29) 


T  =  (T-TJP,  t  =  — ,  y  =  ^. 

tob  h 

In  defining  the  Fourier  number  above,  the  parameter  t*  represents  an  observation  time. 

Equation  (29)  was  treated  numerically  using  a  marching  technique  by  Dr.  Michael 
Khonsari  at  University  of  Pittsburgh.  The  initial  and  boundary  conditions  are  maintained  at 
T  =  0 .  The  time  history  of  the  temperature  distribution  across  the  film  is  presented  in  Figure  28 
for  Br  =  6.75.  As  shown,  temperature  rapidly  rises  in  the  center  of  the  film  at  about  0.06s  and 
becomes  unbounded  when  t  >  0.062s.  The  existence  of  hot  center  temperature  surrounded  by 
cooler  liquid  is  indicative  of  the  onset  of  localization.  This  is  also  apparent  in  the  velocity  profile 
at  t  =  0.062s  shown  in  Figure  28. 


11.4  Optical  Ray  Tracing 

It  seems  reasonable  that  the  last  temperature  profile  of  Figure  28  would  result  in  the 
appearance  of  a  dark  band  at  the  midplane.  However,  a  calculation  is  in  order,  to  show  that  the 
dark  band  accompanies  localization  and  that  the  band  will  not  appear  earlier.  The  previous 
thermal  analysis  yielded  temperatures  at  fifty  equally  spaced  grid  points  within  the  liquid  film.  We 
assume  that  the  temperature  is  uniform  witWn  strata  centered  at  these  points.  Sources  of  light  are 


72 


30  ms  \  t=  60  ms  \  t-  62  ms 


73 


TRANSMISSION  (COUNTS  PER  GRID  POINT) 


placed  at  each  grid  point  across  the  film  at  one  edge  of  the  liquid  film  where  illumination  enters. 
Forty  rays  originate  fi’om  each  source,  equally  distributed  over  an  angular  interval  of  ±0.05  radian 
fi-om  the  centerline  of  the  stratum.  Where  a  ray  intersects  the  boundary  between  strata,  Snell’s 
Law  is  applied.  Contact  with  a  solid  boundary  surface  absorbs  a  ray. 

The  refractive  index,  rio,  of  the  liquid  at  the  test  pressure  and  initial  temperature.  To,  was 
measured  in  a  high-pressure  refi'actometer.  The  Lorenz-Lorentz  equation  can  be  manipulated  to 
yield  an  expression  for  the  refractive  index,  n,  at  other  temperatures  through  the  liquid  density,  p. 
The  Lorenz-Lorentz  equation  [55]  reads 


n^-1 

n^+2 


1 

—  =  constant 

P 


(15) 


The  thermal  expansivity  of  the  liquid  lubricant  was  supplied  by  the  manufacturer  and  corrected  for 
elevated  pressure  from  the  results  of  Yasutomi,  et  al.  [18]. 

The  number  of  rays  exiting  the  film  gap  within  a  stratum  is  considered  a  measure  of  the 
relative  light  transmission  at  the  grid  point  associated  with  that  stratum.  These  values  are  plotted 
in  Figure  28.  It  is  clear  that  the  dark  band  for  which  there  is  relatively  little  light  transmission  is 
not  present  at  60  ms  but  is  fully  developed  at  62  ms  and  although  it  is  a  direct  consequence  of  the 
temperature  profile  the  band  is  associated  with  the  extreme  localization  of  shear  as  indicated  by 
the  velocity  profile  of  Figure  28. 


11.5  Discussion 

Clearly,  the  dark  band  observed  in  the  flow  visualization  experiment  is  due  to  a  thin  hot 
layer  of  liquid  sandwiched  between  two  cooler  layers.  This  hot  layer  represents  a  region  where 
the  shear  deformation  is  localized.  The  greatest  part  of  the  relative  velocity  is  accommodated 
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with  a  small  fraction  (-5%)  of  the  film  thickness.  Notably,  the  localization  occurs  at  shorter  times 
than  observed  in  experiments.  We  believe  that  this  is  due  in  most  part  to  the  inability  of  the 
experimental  technique  to  achieve  true  stress  control  stemming  from  viscous  seal  fiiction  which 
increases  with  sliding  velocity  and  decreases  with  time  of  sliding  or  temperature. 

The  one-dimensional  analysis  presented  here  forces  symmetry  about  the  center  of  the  film. 
The  experimental  geometry  is  finite  in  length  in  the  flow  direction  with  an  inlet  to  the  left  in  Figure 
26.  The  adiabatic  shear  band  near  the  entrance  is  therefore  displaced  toward  the  hotter  stationary 
boundary  at  the  bottom  of  each  print  in  Figure  26  by  the  cooler  liquid  converted  into  the  gap  by 
the  upper  moving  boundary.  This  accounts  for  the  curvature  of  the  band  pictured  in  Figure  26. 

Lubricant  resident  times  in  EHD  are  typically  about  1  ms  or  much  less.  In  this  simulation, 
localization  occurs  at  longer  times  -  on  the  order  of  100  ms.  However,  in  the  dimensionless 
energy  equation  (30),  time  is  scaled  by  h^pCp/Tc.  The  film  thickness,  h,  is  at  least  100  times  thinner 
for  a  lubricated  contact.  Hence,  the  localization  time  would  be  of  the  order  of  10  ps  which  is  in 
the  realm  of  EHD. 

12.  OBSERVATIONS  IN  AN  OPERATING  EHD  CONTACT 

If  light  is  scattered  by  reflection  from  a  mechanical  shear  band,  there  should  be  a  preferred 
scattering  direction  within  the  lubricant  film  which  depends  upon  the  direction  of  shear.  The 
orientation  of  the  second  type  of  band  is  little  changed  by  the  direction  of  shear  and  we  would 
expect  that  the  way  in  which  light  interacts  with  this  type  should  be  little  afiected  by  the  shear 
direction.  In  the  following  discussion,  then,  only  the  first  type  of  mechanical  shear  band  is 
considered. 
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An  optical  EHD  simulator  was  employed  in  this  work  and  is  described  in  Ref  [13],  The 
contact  consisted  of  a  25.4  mm  radius  sapphire  disc  running  over  a  7.6  mm  radius  steel  roller. 
The  roller  had  a  crown  radius  to  ^ve  the  contact  an  ellipticity  of  3.  The  disc  was  loaded  by  an  89 
N  weight  giving  a  contact  size  of  672  x  224  pm  and  an  average  pressure  of  760  MPa.  The 
polyphenyl  ether,  5P4E,  was  used  as  a  lubricant  giving  a  minimum  film  thickness  of  approximately 
1.5pm  at  room  temperature. 

The  original  idea  was  to  illuminate  the  contact  with  a  laser  from  directly  above,  through  a 
piece  of  card  with  a  hole  in  it,  and  observe  the  pattern  reflected  onto  the  underside  of  the  card.  It 
was  hoped  that  light  would  be  reflected  fi'om  the  shear  bands,  expected  to  be  inclined  at  about  19° 
to  the  horizontal,  producing  a  pattern  which  changed  sides  when  the  shearing  direction  was 
reversed.  This  was  not  observed.  Further  examination  of  the  optics  of  the  system  reveals  why.  It 
transpires  that  any  light  which  may  have  been  reflected  firom  the  shear  bands  regardless  of  the 
direction  of  shear  will  strike  the  upper  surface  of  the  disc  at  an  incident  angle  greater  than  the 
critical  one  for  that  interface,  so  the  light  will  be  totally  internally  reflected,  and  could  only  escape 
from  the  disc  edges.  The  refractive  index  of  sapphire  is  1.766.  The  refractive  index  of  the 
lubricant  was  measured  with  a  high-pressure  refract ometer  at  pressures  to  300  MPa  and 
extrapolated  to  760  MPa.  The  value  obtained  is  1.722. 

A  solution  to  this  problem  is  to  illuminate  the  contact  fi'om  an  angle  (Figure  29).  An 
incident  angle  of  about  34°  will  be  refi'acted  so  as  to  strike  a  shear  band  produced  by  negative 
sliding  velocity  approximately  normally,  and  hence  be  reflected  back  into  the  laser.  Consequently 
an  angle  slightly  different  was  used  during  the  experiments.  It  was  hoped  that  this  would  allow 
the  observation  of  light  reflected  back  toward  the  source  in  one  direction  of  shearing  as  illustrated 
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Figure  29.  Schematic  of  light  path  striking  a  shear  band  with  negative  sliding 


in  Figure  29.  Unfortunately,  any  light  that  may  have  been  reflected  was  either  too  difiuse  or 
coincided  with  a  different  brighter  reflection.  However,  it  was  noticed  that  the  intensity  of  light 
escaping  from  the  edge  of  the  disc  (see  Fig.  30)  varied  with  the  direction  of  rotation  and  hence 
shearing.  See  Table  5.  tests  A  through  D.  This  effect  could  be  caused  by  shear  bands,  as 
explained  above  or  may  have  been  due  to  reflections  from  the  inlet  or  outlet  regions  of  the  film,  as 
the  laser  has  a  480pm  diameter  beam  some  light  will  strike  in  front  of  and  behind  the  contact 
patch.  The  inlet  meniscus  is  far  enough  away  to  be  of  no  concern.  However,  the  outlet  region  is 
much  closer  to  the  contact  patch.  The  lubricant  film  cavitates  just  beyond  the  contact  and  most  of 
the  film  is  thought  to  adhere  to  the  surfaces.  If  this  is  the  case  it  will  have  little  effect  on  the 
direction  of  reflected  light. 

To  illustrate  that  this  is  the  case,  the  slide-to-roll  ratio  was  set  negative.  That  is,  rotation 
in  the  same  direction  as  before  results  in  shear  in  the  opposite  one  hence  changing  the  inclination 
of  the  shear  band.  The  higher  intensity  light  was  now  observed  in  the  opposite  direction  of 
rotation  to  before  (i.e.  with  the  outlet  meniscus  on  the  opposite  side  of  the  contact).  See  test  F  - 
Table  5.  The  slide-to-roll  ratio  was  also  set  to  zero  and  no  change  in  intensity  was  observed  when 
the  direction  of  rotation  was  reversed  (test  E).  In  order  to  quantify  this  change  in  intensity  a 
photovoltaic  cell  was  masked,  so  that  only  light  from  the  disc  could  enter,  and  placed  near  the 
disc  edge.  The  output  from  the  photovoltaic  cell  is  only  an  indication  of  the  relative  intensity  for 
rotation  in  opposite  directions  for  each  absolute  value  of  slide-to-roll  ratio,  because  it  proved  to 
be  very  sensitive  to  ambient  lighting  conditions  and  proximity  to  the  disc  edge,  which 
unfortunately  changed  each  time  the  absolute  value  of  slide-to-roll  ratio  was  changed.  The  results 
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Lubricant  Film 


Shear  Band 


Figure  30.  Schematic  of  light  path  striking  a  shear  band  with  positive  sliding 
velocity. 


are  listed  in  Table  5  where  light  intensity  at  the  edge  of  the  sapphire  disc  is  represented  by  the 
photovoltaic  output  in  mV. 

These  results  show  that  the  change  in  observed  intensity  is  not  dependent  on  the  direction 
of  rotation  (i.e.  meniscus  position)  but  on  the  direction  of  shearing.  The  fact  that  the  intensity 
change  can  be  seen  for  both  directions  of  rotation  helps  to  eliminate  the  possibility  that  the 
phenomena  may  be  caused  by  the  geometry  of  the  contact  changing,  or  the  position  and  shape  of 
the  inlet  and  outlet  regions  chan^g  as  the  direction  is  reversed.  While  this  does  not  conclusively 
prove  that  mechanical  shear  bands  exist  in  operating  contacts  (change  in  birefringence  of  sapphire 
with  stress  is  a  possible  effect),  it  does  show  that  there  is  some  mechanism  at  work  within  the 
lubricant  which  is  dependent  on  the  direction  of  shearing  and  effects  the  optical  character  of  the 
lubricant  film. 


80 


Table  5.  Light  Intensities  for  Various  Contact  Positions 


Test 

Rolling 

Velocity 

(m/s) 

Sliding 

Velocity 

(m/s) 

A 

0.228 

0.022 

-0.228 

-0.022 

B 

0.233 

0.013 

-0.233 

-0.013 

C 

0.235 

0.008 

-0.235 

-0.008 

D 

0.237 

0.005 

-0.237 

-0.005 

E 

0.239 

0.000 

-0.239 

0.000 

F 

0.245 

-0.012 

-0.245 

0.012 

Slide-to 
Roll  Ratio 

Intensity  of 
Light  at 
Disc  Edge 
(mV) 

Change  in 
Intensity 
(mV) 

0.096 

51 

0.096 

43 

-8 

0.056 

48 

0.056 

35 

-13 

0.034 

54 

0.034 

37 

-17 

0.021 

40 

0.021 

39 

-1 

0.000 

47 

0.000 

47 

0 

-0.049 

52 

-0.049 

59 

+7 
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Shear  band  analysis  for  lubricants 
based  on  a  viscoelastic  piasticity  modei 

YK  Lee,  J  Ghosh,  S  Bair,  and  WO  Winer 

George  W  Woodruff  School  of  Mechanical  Engineering 
Georgia  Institute  of  Technology,  Atlanta  GA  30332-0405 

Shear  banding  is  an  outstanding  problem  in  lubricant  rheology  where  a  thin  film  of  lubricant 
is  under  high  pressure  and  high  shear  stress.  To  study  such  a  phenomenon,  a  viscoelastic- 
plastic  model  is  proposed.  It  is  postulated  that  shear  bands  appear  when  the  character 
of  the  field  equations  changes  from  elliptic  to  hyperbolic.  The  model  is  derived  based 
on  a  rate  formulation  which  combines  a  MaxweU  fluid  model  and  a  compressible  rate- 
independent  plasticity  model.  The  model  gives  the  constitutive  relation  of  a  compressible 
viscoelastic-plastic  fluid  for  which  the  assumption  of  Stokes  condition  becomes  unnecessary. 
It  also  accounts  for  the  elastic  coupling  effect  of  both  the  viscous  and  the  rate-independent 
behavior  of  a  lubricant.  It  provides  a  novel  feature  that  the  development  of  shear  bands 
can  be  tracked  rather  than  being  determined  after  they  are  fully  developed.  Hence,  the 
process  of  shear  banding  may  be  articulated.  The  analysis  focuses  on  identifying  necessary 
and  sufficient  conditions  which,  if  achieved  in  a  deformation,  would  produce  a  change  in 
character  of  the  governing  field  equations.  Through  the  use  of  Drucker-Prager  criterion 
as  the  component  that  governs  the  rate-independent  behavior  of  a  lubricant,  the  analysis 
gives  results  that  are  in  accord  with  experimental  observations.  Parametric  results  are 
given  graphically. 
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INTRODUCTION 

The  search  for  a  general  rheologiceJ  constitutive  equation 
for  liquid  lubricant  at  high  pressures  has  been  a  research 
focus  for  some  time.  It  has  been  observed  that  under  a 


high  pressure  and  shear  stress,  a  lubricant  may  exhibit 
a  nonlinear  stress  strain-rate  behavior,  whereas  the  same 
lubricant  usually  follows  a  linear  Newtonian  response  be¬ 
low  a  critical  shear  stress  (Bair  smd  Winer,  1992).  The 
deviation  from  the  linear  response  marks  a  character¬ 
istic  change  of  the  lubricant.  It  begins  to  behave  like 
a  viscous  solid  because  there  exists  a  limiting  value  of 
sheeir  stress  at  which  the  response  becomes  strain-rate 
independent.  Bair  and  Winer  (1992)  suggested  a  nonlin¬ 
ear  Maxwell  model,  with  shear  modulus  taken  into  the 
time  derivative,  which  included  the  effect  of  pressure  on 
viscous-plastic  transition.  In  a  more  recent  work,  Bair 
et  al.(1993)  reported  the  observation  of  short-lived  shear 
bands  in  a  lubricant.  It  was  observed  that  the  first  shear 
bands  appear  at  the  onset  of  nonlinear  behavior  marked 
on  a  typical  stress  strain-rate  curve.  An  increase  in  the 
8he2ir  stress  above  this  point  is  accompanied  with  a  prolif¬ 
eration  of  shear  bands.  The  inclination  of  the  shear  bands 
appears  to  be  independent  of  the  pressure  and  the  shear 
rate.  Bair  and  Winer  (1993)  also  showed  that  shear  b2md 
angles  can  be  predicted  by  using  the  Mohr-Coulomb  crite¬ 
rion.  The  orientation  of  shear  bands  after  the  localization 
has  occurred  could  be  obtained  from  the  yield  trajectory 
by  using  a  proper  planer  stress  transformation.  However, 
there  still  exists  a  need  for  conditions  that  could  predict 
the  onset  of  shear  banding. 
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In  this  study,  a  viscoelastic  plasticity  model  is  derived 
based  on  a  rate  formulation  which  combines  a  Maxwell 
fluid  model  rmd  a  compressible  plasticity  model.  The 
model  accounts  for  the  elastic  coupling  effect  of  both  the 
viscous  and  the  rate-independent  behavior  of  lubricants. 
With  this  model,  the  necessary  and  sufficient  condition 
for  the  onset  of  shear  banding  is  presented  following  the 
standard  procedure  developed  by  Lee  (1989,  1992),  in 
which  shear  banding  is  regarded  as  a  geometric  instabil¬ 
ity,  and  the  onset  of  such  an  instability  is  defined  at  the 
instant  when  the  character  of  the  governing  field  equa¬ 
tions  ceases  to  be  elliptic. 

In  what  follows,  we  first  provide  a  framework  for  the 
rate  formulation  of  viscoelastic  plasticity.  This  framework 
can  be  used  to  simulate  the  combined  rate-dependent 
and  rate-independent  deformation  behavior  in  either  liq¬ 
uid  or  solid  polymers  while  analysis  can  adopt  various 
approaches  to  plastic  dissipation  potential,  we  choose 
the  Drucker-Prager  (1952)  criterion  to  account  for  the 
pressure-dependent  deformation  behavior  that  a  lubricant 
may  exhibit.  This  criterion  is  chosen  because  it  is  math¬ 
ematically  simple  and  yet  physically  plausible  for  liquid 
medium  which  cannot  sustain  any  tension.  It  must  be 
emphasized  that  the  choice  of  such  a  plastic  dissipation 
potential  does  not  restrict  in  any  manner  the  develop¬ 
ment  of  the  procedure  for  shear  band  analysis.  It  simply 
provides  a  means  to  include  the  effects  of  hydrostatic  pres¬ 
sure  that  is  important  to  the  problem  at  hand.  Of  course, 
the  ultimate  justification  of  using  a  specific  plastic  dissi¬ 
pation  potential  and  its  domain  of  validity  depend  on  the 
ability  of  the  resulting  model  to  predict  experimental  ob¬ 
servations. 

VISCOELASTIC  PLASTICITY  MODEL 
Additive  decomposition 

One  of  the  common  assumptions  in  various  inelasticity 
theories  is  the  decomposition  of  strain  into  an  elastic  and 
an  inelastic  pturt,  although  details  of  decomposition  may 
vary  from  one  material  to  another,  depending  on  the  be¬ 
havior  that  needs  to  be  emphasized.  To  describe  the 
flow  behavior  of  lubricants  and  to  account  for  both  rate- 
dependent  and  rate-independent  behavior  at  high  pres¬ 
sure,  we  assume  that  the  total  rate  of  deformation  ten¬ 
sor  iij  =  {vij  -b  Vj,i)/2  for  a  given  velocity  field  Vi  can 
be  decomposed  into  a  viscoelastic  component  iff  and  a 
rate-independent  plastic  component 

=  +  (1) 

Note  that  we  will  use  ” ,”  before  a  subscript  to  denote  pair- 
tial  derivatives  throughout  the  paper.  In  what  follows,  we 
first  derive  an  expression  for  eJ)*  based  on  thermodynamic 
potentials  and  then  for  based  on  a  rate-independent 
dissipation  potential  using  the  norm2dity  flow  rule. 


Viscoelastic  formulation 

Various  models  could  be  used  to  model  the  response  of 
viscoelastic  liquids  (Lemaitre  and  Chaboche,  1990).  In 
this  investigation,  we  adopt  a  Maxwell  fluid  model  for 
its  simplicity  and  yet  its  accountability  in  explaining  the 
e^qperimenteilly  observable  relaxation  of  stress  (Bair  and 
Winer,  1993).  It  is  assumed  that  the  total  viscoelastic 
strain  is  the  sum  of  an  elastic  and  anelastic  part,  so  that 


f?'.*  =  £?.  -I- 

~  ‘■tj  • 

The  elastic  component  of  the  strain  is  given  by 


£ii  =  P 


dip* 

dau 


and  the  anelastic  component  in  a  rate  form  by 

dtp* 


£?r»  = 


dcTij' 


(2) 

(3) 

(4) 


where  ip*  and  tp*  are  the  thermodynamic  potential  and 
the  dual  dissipation  potential,  respectively.  Both  poten¬ 
tials  are  functions  of  the  stress  state. 


For  a  linear  theory,  the  S2une  thermodynamic  potential 
can  be  used  as  for  the  isotropic  linear  elasticity.  We  thus 
write 

where  p  is  the  mass  density  of  the  material,  v  the  Pois¬ 
son’s  ratio,  E  the  elastic  modulus,  and  ?  the  stress  tensor. 
SimilMly,  the  dud  dissipation  potential  tp'  is  defined  by 
a  positive  definite  quadratic  form  which  includes  two  new 
coefficients  ri  and  T2  and  is  given  by 


The  coefficients  n  and  r2  are  characteristic  of  the  vis¬ 
cosity  emd  they  may  be  determined  by  shear  and  uniax¬ 
ial  compression  tests.  The  viscoelastic  rate  of  strain  can 
then  be  determined  by  combining  Eqs  (2)-(6).  The  result 
is  given  by 


where  Sij  is  the  Kronecker  delta.  The  expression  in  the 
first  bracket  is  the  elastic  component.  The  one  in  the 
second  bracket  takes  into  account  the  effects  of  viscosity 
which  result  in  a  time  dependent  response  of  a  viscoelastic 
medium.  Hence  Eq  (7)  can  be  cast  into  the  standard  form 


iij  =  SijkiCTki  +  Rijkio-ki, 


(8) 
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where  Sijti  is  the  elastic  compliance  tensor  and  Rijki  the 
viscous  creeping  tensor.  Inverting  Eq  (8),  and  then  ex¬ 
pressing  the  results  in  terms  of  Lame’s  constants  A  and 
fi,  we  have 

Cij  =  [2pe?/  +  A«yej;|] 

or,  in  short 

~  ■AijklO'kl,  (10) 

where  the  component  of  the  viscoelastic  modulus 

^ijkt  =  K^ikSjl  -I-  SilSjk)  -h  X(8ij6kt)  (11) 

and  the  components  of  the  viscous  relaxation  tensor 


■^jkl  —  p  {Sik^jl  +  ^il^jk)  +  A*(iyitf)  (12) 


with 


E 

2(1  +  1/) 


and 


A  = 


Ev 

(l  +  i/)(l-2i/)' 


A'  = 


l-2i/ 


Ln 


(13) 

(14) 


It  can  be  seen  that  the  expression  in  first  bracket  of  Eq 
(9)  is  of  a  form  identical  to  the  generalized  Hooke’s  law 
for  an  isotropic  medium,  except  that  the  components  are 
of  the  viscoelastic  rate  of  deformation  tensor  instead  of 
being  the  elastic  components  of  the  total  deformation  rate 
tensor.  The  expression  in  the  second  bracket  accounts  for 
the  time-dependent  viscous  effects  such  as  relaxation  and 
creep  that  depends  on  the  instantaneous  state  of  stress 
in  the  material.  The  physical  interpretation  of  Eq  (9) 
is  as  follows;  The  increment  of  the  stress  at  any  instant 
of  time  depends  on  an  instantaneous  Eookean  response 
of  the  viscoelastic  straun  increment.  The  instantaneous 
stress  increment  includes  all  rate  and  time  dependent  ef¬ 
fects,  and  during  the  increment,  the  viscous  relaxation 
and  other  time  dependent  behavior  of  materiad  occurred 
in  the  previous  stress  state  is  removed. 


Creep  and  relaxation  responses 


The  application  of  Eq  (7)  to  a  plane  shear  test  gives 


1  + .  1  + 


(15) 


In  response  to  relaxation  test  ejl  =  0,  the  integration  of 
Eq  (15)  leads  to 


<^12  =  ^12  ®CP 


(16) 


where  instantaneous  shear  stress  at  time  f  =  0. 

Thus  Ti  can  be  physically  interpreted  as  the  relaxation 


time  in  a  plane  shear  test  and  can  be  determined  from 
Eq  (16). 


•  As  indicated  by  Eq  (1),  the  totad  rate  of  deformation  is 
the  sum  of  the  viscoelastic  part  and  the  rate-independent 
part.  The  volumetric  creep  and  relaxation  behavior  of  a 
lubricant  under  a  hydrostatic  compression  can  be  studied 
by  exaunining  the  following  rate  equation 


l-2i/.  ^ 


fl  +  v 
\En 


Vkk- 


(17) 


At  t  =  0, 

£2U0)  =  i^^r«(0)  (18) 

is  the  initial  condition,  which  expresses  an  instantaneous 
elasticity.  Through  some  algebraic  manipulation,  the 
creep  function  under  a  step  load  function  akk{t)  =  1(f) 
can  be  found  by  solving  Eq  (17),  amd  the  result  is  given 
by 

4t(f)  =  A  +  [Afi(f)  +  B]f  (19) 

where  A  =  (1  —  2j/)/E,  B  =  (1  +  v)/Eti  —  ZvjEr^^  and 
£(f )  indicates  the  Dirac-delta  function,  which  is  defined  as 
a  function  with  singularity  at  t  =  0  and  6(t)  =  0  for  f  <  0, 
and  f  >  0.  The  relaxation  behavior  is  described  by  the 
relaxation  function  o-*t(f)  which  renders  the  deformation 
prescribed  by  e^Kt)  =  1(f).  This  relaxation  function  can 
be  found  by  solving  Eq  (17)  again.  The  result  gives 

<^kk{t)  =  jexp  +  ;^^(0  (20) 


The  quantity  —A/B  can  be  interpreted  as  the  relaxation 
time,  a  characteristic  of  stress  relaxation.  If  the  relax¬ 
ation  parameter  ti  is  known,  the  parameter  can  be 
determined  using  the  above  equation. 

Note  that  the  creep  function  and  the  relaxation  func¬ 
tion  are  characteristic  only  of  the  viscoelastic  properties 
of  the  lubricant  under  test.  These  two  properties  are  de¬ 
fined  in  terms  of  the  deviation  of  pressure  from  a  reference 
pressure. 

General  viscoelastic  plastic  formulation 

Noting  the  relation  in  Eq  (1),  We  rewrite  Eq  (10)  as 


d’y  —  Gyn(CiJ  —  ^j)  —  Aijkl<Tkl-  (21) 

This  equation  represents  a  general  viscoelastic  plastic 
flow  relation  with  a  coupled  rate-dependent  and  rate- 
independent  behavior.  It  is  clear  that  at  any  instant 
of  time  we  can  calculate  the  stress  rate  for  a  given  rate 
of  deformation  and  stress  state,  provided  that  the  rate- 
independent  plastic  component  of  the  strain  rate  is 
known.  This  cem  be  determined  once  the  form  of  a  suit¬ 
able  rate-independent  plastic  potential  ^  (iTy )  is  adopted. 
The  derivation  can  follow  that  in  a  rate-independent  plas¬ 
ticity  theory. 
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Loading  function  and  flow  rule 

Assume  there  exists  a  loading  function  given  by 

=  (22) 

Equation  (22)  represents  a  surface  in  the  stress  space. 
Alternately,  we  write 

-  X{W^)  =  0.  (23) 


Letting  4>ij  =  d4>ld<Tij  zmd  replacing  A  in  Eq  (26)  from 
Eq  (28),  and  then  proceeding  to  Eq  (21),  we  obtaiin 

^ij  =  Cijiiiitl  —  (Caie4®e4  ~  Aated(^cd) 

-  Aijkicrku  (30) 

or,  in  a  more  familiar  form 


where  is  the  effective  stress  state  and  the  response 
to  the  effective  stress  is  given  by  x(W^)  which  depends 
on  the  plastic  strain  energy  density  W^,  a  function  of 
accumulated  plastic  strain  ef^- .  The  rate  change  of  plastic 
strain  energy  is  given  by 

VH’  -  =  <re,e5,,  (no  sum  oneg)  (24) 

where  a-e^,  are  the  equivalent  stress  and  the  equivalent 
plastic  strain  rate,  respectively.  During  deformation,  the 
consistency  condition  (Prager,  1949)  which  simply  states 
that  the  loading  state  must  trace  its  own  stress  strain 
curve  requires  that 

F  =  ^-X  =  0.  (25) 

Instead  of  treating  the  normality  rule  as  an  assumption, 
Lee  (1993)  showed,  through  the  use  of  Legendre  dual 
transformation,  that  Drucker’s  hypothesis  (1964)  of  path 
independence  in  the  small  is  a  characteristic  of  the  addi¬ 
tive  decomposition  assumption  of  strain  rates.  It  was  also 
shown  that  normality  of  plastic  strain  rate  to  the  loading 
function  is  a  consequence  of  the  path  independence  in  the 
small  even  when  a  material  instability  occurs. 

The  normality  flow  rule  associated  with  the  loading 
function  F((t,j,  W^)  is  given  by 


^ 


where  A  is  a  parameter  expressing  the  magnitude  of  the 
plastic  component  of  the  rate  of  deformation  and  can  be 
determined  from  the  consistency  condition  Eq  (25).  Sub¬ 
stituting  Eq  (21),  together  with  the  normality  flow  rule, 
into  the  consistency  condition,  we  have 


dx  dwr  i  .  „ 

dwp  d^ij  dcij 
Solving  for  A,  we  obtain 
I  d<i>  , 

A  =  f  iPijkl^kl  -  AijkKTkl)  , 

where  the  scalar 

d4>  d<t>  dx  dW»  d<l> 

~  daij'^^'dcTki  dWP  de^ij  dffij- 


(27) 

(28) 

(29) 


&ij  —  Cijkt^kl  “  Vijkl^^kh 

(31) 

where 

Cijn  =  Cijij  —  Cijti 

(32) 

and 

Vijjci  =  Aijti  — 

(33) 

with 

Clj^i  —  ~^<l>ai^ed^ai^hj^ek^dU 

(34) 

Aij^i  —  p  ^ab^edSai^hj^ek^dlt 

(35) 

^ij  =  4>ij  + 

(36) 

A' 

’ 

(37) 

^ijkt 

and  Aijki  are  defined  before. 

FIG.  1  Stress-str2dn  curve  for  a  viscoelastic  plastic  medium. 

The  interpretation  for  Eq  (31)  is  shown  in  Fig.  1  graph¬ 
ically.  The  coefficient  Cfj^j  can  be  viewed  as  the  com¬ 
ponents  of  a  corrected  elastic  modulus  which  results  in 
the  stress  increment  taken  up  for  rate  independent  plas¬ 
tic  strain,  when  applied  directly  to  the  tottd  rate  change 
of  deformation.  Similairly,  can  be  interpreted  as 
the  components  of  a  corrected  time-dependent  relaxation 
modulus  which,  when  applied  to  the  stress  state  at  the 
instant  the  increment  of  deformation  is  applied,  results 
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in  rate-independent  plastic  strain  that  cannot  relax  with 
time.  Therefore,  we  can  define  a  viscoelastoplastic  mod¬ 
ulus  Cijn  given  by  Eq  (32)  which  provides  an  apparent 
stress  rate  dependent  on  a  Hookean  response.  This  ap¬ 
parent  stress  rate  has  to  be  corrected  for  effective  viscous 
relaxation  and  creep,  by  subtracting  the  stress  rate  result¬ 
ing  from  the  above  mentioned  effects,  dependent  on  the 
stress  state  at  the  instant  the  increment  of  deformation  is 
applied.  Thus  Vijn  in  Eq  (33)  is  an  effective  viscoelastic 
relaxation  tensor  which  applies  to  the  current  stress  state 
and  provides  the  effective  time  dependent  effects. 

Symmetry  property 

The  above  formulation  is  based  on  the  additive  decom¬ 
position  of  the  rate  of  deformation  tensor  as  indicated  in 
Eq  (1).  Such  an  assumption  leads  to  a  quasilinear  rela¬ 
tionship  between  Cij  and  iij.  Generally  speaking,  this 
linearity  has  been  foimd  to  be  true  for  thermodynamic 
systems  sufficiently  close  to  equilibrium.  Of  course  the 
ultimate  justification  of  this  property  relies  on  the  abil¬ 
ity  of  the  representing  equations  to  predict  experimental 
observations.  For  deformation  to  be  path  independent 
in  the  small,  Cyn  must  be  symmetric  (Lee,  1993).  This 
symmetry  property  always  exists  because  both  ^ijkl  and 
Cljj.1  are  symmetric.  To  have  a  symmetric  stress-strmn 
response,  however,  we  also  require  that  the  coefficient  of 
stress  in  Eq  (31)  be  symmetric  as  well.  This  can  be  made 
by  imposing  a  symmetric  condition  on  Thus,  we 

require  that 

^ah^cdiaiibjSekSdl  =  iai^ci^akhl^dSdj  ■  (38) 

Substituting  ^ab  and  4ed  from  Eqs  (36)  and  (37),  we  ob¬ 
tain 

/A  X'  \ 

“  V  j  “  <i>uSki<i>i}]  =  0.  (39) 

Since  the  terms  in  the  bracket  cannot  be  zero  should  any 
deformation  occur,  we  can  only  require  that 


Substituting  from  Eqs  (13)  and  (14),  we  conclude  that 
either  n  =  0  or  1/t2  =  0. 

We  recall  that  ti  represents  the  characteristic  relax¬ 
ation  time  for  a  shear  deformation;  it  cannot  vanish  be¬ 
cause  lubricants  usually  exhibit  a  marked  rate-dependent 
behavior  in  shear.  Thus  the  only  logical  choice  is  that 
1/t2  must  vanish.  As  a  result  of  this  condition  and  the 
discussion  on  the  relaxation  constants  in  a  previous  sec¬ 
tion,  the  complete  relaxation  behavior  of  lubricant  is  de¬ 
scribed  by  a  single  relaxation  par2Uiieter  Tj.  Note,  how¬ 
ever,  that  we  will  continue  to  have  different  characteristic 
time  for  the  two  deformation  modes,  namely  shear  defor¬ 
mation  and  hydrostatic  compression.  From  the  symmetry 
requirements,  we  find  from  Eqs  (36)  and  (37)  that 
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Physical  interpretation  of  F 

Alternatively,  we  can  solve  A  by  combining  Eqs  (24),  (25) 
and  (26);  the  result  is 

k  =  (42) 

where  the  inverse 

Since  is  a  scalar  measure  of  the  effect  of  the  stress 
tensor  on  the  material,  we  infer  that  4>{pi})  is  a  flow  sur¬ 
face  which  distinguishes  the  nonlinear  rate-independent 
viscous  behavior  (which  is  termed  plastic  before)  from 
the  rate-dependent  one.  We  can  represent  this  surface  by 
an  equivalent  stress,  <r,q.  Correspondingly,  there  is  some 
equivalent  plastic  strain,  whose  increment  produces  a  rate 
change  in  the  rate-independent  plastic  strain  energy  asso¬ 
ciated  with  the  material.  Hence,  from  Eq  (24)  it  is  clear 
that 

=  CijSHij  (44) 

Noting  that  ^  =  Ctq,  we  write 

~  (45) 

where  the  instantaneous  plastic  tangential  modulus, 

2^§j  =  dej/cfj,  represents  the  slope  aJongthe  equivalent- 
stress  equivalent-plastic  strain  curve  of  a  material.  This 
modulus  represents  the  hardening  character  of  the  lubri¬ 
cant.  Using  Eqs  (11),  (36)  and  (43),  we  obtain  from  Eq 
(29)  the  expression  for  F, 

r  =  ^  +  (46) 

Similsir  to  S  physicaJly,  F  represents  the  effect  of  the  mate¬ 

rial  specific  hardening  response  to  the  instantaneous  state 
of  stress  in  the  material. 

With  the  above  results  and  through  some  manipula¬ 
tions,  one  can  show  that  coefficients  Cijki  and  Vijki  in  Eq 
(31)  reduce  to 

Cijki  =  (d{hik6ji  -I-  SitSjk)  +  XSjjSki 

_  ifl^ak^cdSaiSbjSckhdl 
l/0*  +  4>mn^rt6mr6n3 
and 

Vijki  =  ld'{hikSji  +  6ii6jk)  +  X'SijSki 

_  ^n'^ai^cdSaiSjjSckSdl 
l/^*  ^mn^rihmrhns 


4>ij  =  4>ii  ■ 


(41) 


(48) 
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respectively,  where  =  2fiP. 


In  the  past  there  has  been  an  inconclusive  discussion 
relating  to  the  validity  of  Stokes  condition,  and  its  out¬ 
come  on  the  modeling  behavior  of  Newtonian  fluids.  It  is 
clear,  from  the  viscoelastic  plastic  formulation,  that  there 
is  no  need  to  invoke  the  Stokes  condition  in  order  to  ded 
with  the  constitutive  modeling  of  lubricants. 

Example  for  die  flow  potential 

So  far  we  have  proposed  a  general  formulation  for  the 
viscoelastic  plasticity  of  lubricants  without  having  dis¬ 
cussed  a  specific  form  for  the  rate-independent  plastic 
flow  potential.  Bowden  and  Jukes  (1972)  reported  the 
yield  behavior  of  lubricant  and  liquid  polymers;  their  re¬ 
sults  indicated  that  yielding  is  pressure  dependent.  They 
suggested  the  use  of  a  modified  von  Mises  criterion  to 
describe  such  a  behavior  for  polymers.  Their  criterion, 
in  fact,  is  termed  as  the  Drucker-Prager  criterion  (1952). 
It  is  often  used  in  soil  mechanics,  for  which  soil  is  usu¬ 
ally  treated  as  a  medium  that  cannot  support  a  tension 
force,  a  characteristic  exhibited  by  all  liquids  as  well.  In 
our  analysis,  we  have  adopted  the  same  criterion,  but  we 
write  it  in  an  alternative  form  in  terms  of  a  mean  pressure 
P 

=  VJi  +  aP,  (49) 

where  is  the  second  invariant  of  the  stress  deviator 
tensor  Sy,  and  P  is  positive  for  compression.  For  a  rigid 
incompressible  plasticity  model,  a  is  boimded  by  ±-\/Z/2. 
When  a  is  zero  ,  Eq  (49)  reduces  to  the  von  Mises  crite¬ 
rion. 


Clearly,  the  Drucker-Prager  criterion  provides  for  a  lin¬ 
ear  increase  in  the  shear  yield  stress  with  the  mean  hy¬ 
drostatic  pressure  P.  This  type  of  linear  variation  of  the 
yield  stress  with  pressure,  has  been  found  to  be  obeyed 
by  lunorphous  glassy  polymers  (Bowden,  1973).  The  di¬ 
mensionless  constamt  of  proportionality  a  is  a  material 
property.  Its  value  has  to  be  determined  by  fitting  exper¬ 
imental  results.  It  csin  be  shown  from  Eq  (49)  that 


4>ii 


(50) 


EQUILIBRIUM  REQUIREMENT 

A  suitable  freune  invariant  stress  rate  must  be  used  with 
the  constitutive  relation  presented  in  the  previous  sec¬ 
tions.  The  selection  of  a  frame-indifference  stress  rate  is 
by  no  means  conclusive.  To  provide  a  &ame-indifference 
flow  rule  we  replace  the  stress  rate  dy  in  Eq  (31)  by  the 
Jaumaim  stress  rate  dy.  The  relation  between  these  two 
stress  rates  is  given  by 

^y  =  d’y  (Titn^mj  ~  (51) 

where  tiy  is  the  spin  tensor  and  Wij  =  (u.  j  —  t;j,i)/2. 


Since  deformation  is  considered  as  a  rate  problem,  the 
requirements  of  mechanical  equilibrium  must  be  applied 
not  only  to  the  instantaneous  state  of  stress,  but  zdso  to 
their  rate  change.  The  derivation  of  stress  rate  equations 
of  equilibrium  can  be  found  elsewhere,  e.g.,  see  the  work 
by  Lee  (1988).  The  requirement  states  that 

~  =  0.  (52) 

This  equation  gives  the  complete  description  for  the  equi¬ 
librium  state  of  the  material.  It  governs  the  stress  field  ir¬ 
respective  of  the  magnitude  of  deformation,  and  material 
structure.  Satisfaction  of  equation  (52)  not  only  implies 
total  stress  equilibrium,  but  edso  assures  that,  given  an 
equilibriated  stress  field,  equilibrium  is  maintained  in  the 
presence  of  time  dependent  loading.  Note  that  the  con¬ 
vective  rate  of  change  of  stress  (rkijVj,k  in  Eq  (52)  van¬ 
ishes  only  in  the  case  of  homogeneous  deformation,  and 
thus  &ijj  alone  may  not  represent  the  material  derivative 
of  the  total  stress  equilibrium. 

SHEAR  BAND  ANALYSTS 
Velocity  field  equations 

The  method  of  characteristics  has  proven  to  be  a  fruit¬ 
ful  means  of  representing  the  sheaur  bamds  in  plastic  flow 
problems.  Since  the  characteristics  of  the  pairtial  dif¬ 
ferential  equation  of  hyperbolic  t)fpe  may  indicate  paths 
upon  which  discontinuities  will  propagate,  it  is  logical  to 
assume  that  these  characteristics  represent  shear  bands. 
This  concept  has  been  used  by  Prandtl,  Nadai,  and  other 
investigators  in  studying  slip  lines  in  plastically  deformed 
bodies.  The  approach  to  be  presented  here  is  a  standard 
procedure  for  shear  band  analysis,  at  least  in  the  context 
of  elastoplasticity  (Lee,  1989,  1992). 

In  what  follows,  we  consider  a  plane-strain  deformation 
in  rectangular  cartesian  coordinates,  and  the  plane  nor¬ 
mal  to  the  z-axis,  so  that  the  independent  space  variables 
are  x,  and  y,  df  dz  being  a  null  operator.  The  non-zero 
deformation  rate  and  spin  components  are  Sxzi  £yy>  Cxy> 
and  tixy.  Thus,  in  terms  of  velocity  components,  the 
equilibrium  equations  (52)  reduce  to 

+  ^ry,y  ~  ~  <7xx,yVy,r 

~^xy,xVx,y  “  ^xy,yVy,y  =  0  (53) 

<7xy,x  +  (Zyy^y  —  0’xx,xVx,x  ~  ^xy.yVy.x 

— (Tyy^xUx.y  ~  ®'yy,yVy,y  —  0.  (54) 

The  constitutive  relationship  ciin  be  written  in  a  matrix 
form  for  plane  strain  deformation  as 


<7xx 

■  hn 

612 

613 

/  1 

621 

622 

623 

'  A 

ffxy 

bsi 
.  ^41 

^32 

642 

CO  CO 
CO  ^ 
•0  *0 

\  *'y.s 

\  (Vx,y  +  t'y,x)  J 
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an 

012 

Oi3 

014  ' 

(  0** 

021 

O22 

O33 

024 

J  <^yy 

031 

O32 

O33 

O34 

I  <^*y 

041 

O42 

043 

O44  . 

[  ffxx 

Denoting  the  shear  stress  by  T,  i.e.,  and  the  first 

normal  stress  difference  by  2JV,  i.e.,  N  =  (<r**  —  ayy')f’l, 
and  then  combining  the  equilibrium  equation  (53),  the 
constitutive  equation  (55)  in  terms  of  the  Cauchy  stress 
rate  by  using  Eq  (51),  we  obtain 

^iif*,**  +  [^13  +  ^31  +  T]vx,tv  +  [633  — 

+[613  —  T^Vy^gg  +  [612  +  ^33  + 

+[^11,*  +  ^31,2  — 

+[^12,*  +  632,2  ~ 

+[613,*  +  633,5  ~  (^*1,5  “  <^»y,»)/2](t'*.y  +  «y,r) 

— [ail,r  +  a31,y]o’iz  —  Wl2,x  +  a32,y]o'yy 
— [ai3,i  +  a33,y]0’zy  “  [ai4,z  +  a34,y]<7*, 
“flllO"**,*  ~  <*31^^**,y  fll2^yy,z  ~  ®32®’yy,y 
— <ll30'zy,z  ~  033^^zy,y  ~  OU^Zi.z  ~  fl34®'zz,y 

=  0.  (56) 

Similarly,  the  combination  of  Eqs  (51),  (54)  and  (55)  gives 


Defining  new  variables  u,-,  such  that 

til  —  t^Z.Z,  tl2  —  Vx,y,  U3  —  Vy^x,  tl4  —  Vy^y.  (58) 

Equations  (56)  and  (57)  may  now  be  expressed  in  terms 
of  the  new  variables  xn,  and  the  result  can  be  cast  into  a 
canonical  form 

(i,i  =  l,2,3,4).  (59) 

We  may  posit  that  there  are  linesur  combinations  of  Eq 
(59)  such  that  the  result  is  a  common  directional  deriva¬ 
tive  of  ti{.  That  is,  the  expression 

=  (60) 

is  directional,  where  Vi  is  the  direction  of  the  direc¬ 

tional  derivative.  This  derivative  is  along  a  set  of  curves 
i9(x,y,  u,-)  =  constant,  so  that  the  directions  are  normal 
to  the  gradient  of  Thus  we  write 

ViPij^,x  +  =  0.  (61) 

Denoting  the  slope  of  the  characteristic  curves  by  A,  thus 
A  =  (62) 

and  Eq  (61)  becomes 


63lt'z,zz  +  [633  +  621  —  -I-  [623  —  T]t)x,yy 

+[633  +  ./V^)]Vy  ,XX  •f  [632  +  623  +  T\Vy^gy  +  l>22X)y,yy 

+[631, z  +  621,5  ~ 

+[632.z  +  622,5  “  ^■yy.yM^x.y 

+[633,z  +  623,5  -  ^(<^yy.*  -  <^zy,y)](t'z.y  +  Vy,g) 

—[031, z  +  a21,5]o’zz  —  [a32,z  +  a22,y]o'yy 
~[a33,z  +  023,y]4^zy  —  [a34,z  +  a24,5]<^iz 
— fl31<r*z,z  ~  fl21®’zz,5  ■  ■  fl32®’55,z  ~  022^55,5 
— a33<Tj,y,s  —  a230’zy,y  ~  034®'ii,z  ~  a24<^zz,5 

=  0.  (57) 

Equations  (56)  and  (57)  are  the  velocity  equilibrium  equa¬ 
tions  of  the  deformation,  and  govern  the  instantaneous 
spatial  dependence  of  the  velocity  field.  Thereby  they 
provide  a  quasi-linear  model  for  the  entire  deformation 
process. 

Characteristic  analysis 

We  now  examine  the  conditions  under  which  the  velocity 
equilibrium  equations  (56)  and  (57)  may  be  integrated  on 
characteristic  surfaces. 


[g.-;-AP.;]i/.  =  0.  (63) 

For  Vi  to  have  a  non-trivial  solution,  the  determinant  of 
the  coefficients  of  Vi  must  vanish.  Hence, 

det|gy-APy|  =  0.  (64) 

Noting  that  coefficients  a,j,  and  6y  are  symmetric  we 
deduce  from  Eq  (64)  a  quartic  equation  in  A  and  the 
result  is 

A^  [613(613  —  T)  —  611(633  -b  A^)] 

-fA®  [2(611623  -  612613) 

+(611  ■+  612  -b  2b33)T  -b  2613N’] 

-bA^  [612(612  —  2633)  —  611622  —  2613623  —  3(613  +  623)7’] 

-b  A  [2(6136  22  —  612623) 

+(622  +  612  +  2633)r  —  2623iV'] 

+  [623(632  —  7’)  —  622(633  —  N)] 

=  0.  (65) 

Denoting  the  principal  components  of  the  quantities 
Cij,<t>ij,  and  ^ij  by  a  single  subscript,  i.e.,  <ri,..,  ^1,.., 
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etc.,  we  obtain  the  coeflScients  bij  in  terms  of  principal 
coordinates;  the  result  gives 

f>33  =  fi 

=  (i-t)r- 

6i3  =  0 

623  =  0 

643  =  0  (66) 


(72) 

This  gives  the  critical  value  of  N  at  which  shear  band¬ 
ing  becomes  possible.  However,  one  should  note  that  this 
condition  alone  is  not  sufficient.  The  inequality  repre¬ 
sented  by  the  second  condition  in  Eq  (70)  must  also  hold 
for  shear  bamds  to  occur.  Alternatively,  this  condition 
can  be  checked  by  calculating  A  directly  because  coeffi¬ 
cient  a  is  always  positive  for  a  nonzero  stress  state,  and 
the  values  of  A  are  real  if  6  <  0. 

For  a  given  value  of  N,  the  corresponding  critical  value 
off*  can  be  found  by  solving  Eq  (72).  The  result  gives 


where 

r*  =  — -b  4^1^! -f- ^2^2  +  ^3^3-  (67) 

Expressing  in  terms  of  principal  quantities  and  replacing 
N  hy  N  which  is  defined  as  N  =  N/ft,  we  recast  Eq  (65) 
in  the  form 

a  (l-b  JV)  A** -t- 26A*  -b c (l  -  iV)  =  0  (68) 

with 


a  =  (l_,/)r-(l-2i/)^? 

6  =  (1  -  i/)r*  -b  ^1^2  ~  (1  ~  +  ^1) 

c  =  (1  -  i/)r*  -  (1  -  2i/)^i.  (69) 

For  JV  <  1,  the  necessary  and  sufficient  condition  for  the 
onset  of  shear  banding  is  that 

(i)  |JV|  =  [(ac-6V«c]^/=' 

and 

(ii)  6  <  0  (70) 

hold  simultaneously.  The  slopes,  i.e.,  the  chuacteristic 
directions  of  shear  bands,  can  then  be  found  by  solving 
Eq  (68).  The  result  gives 


A  =  ± 


(71) 


RESULTS  AND  DISCUSSION 


Critical  values  of  N,  F*,  and 

Inserting  Eq  (69)  into  Eq  (70)  leads  to 

1-^1  =  I  (^1  -  ^2)  ~  (^1  +  ^2)  + 


Corresponding  to  FJ^,  we  can  determine  the  critical  value 
of  from  Eqs  (46)  for  an  instantaneous  stress  state. 
Corresponding  to  /?*,.,  there  exists  a  critical  plastic  tan¬ 
gential  modulus  of  the  lubricant,  below  which 

sheeir  bamds  may  not  occur.  Since  Per  is  also  a  function 
of  stresses  which,  in  turn,  will  determine  the  inequality  in 
Eq  (70).  Thus,  the  critical  plastic  tangential  modulus  for 
the  onset  of  shear  banding  is  not  a  constemt.  It  is  stress 
and  thus  problem  dependent. 

The  criticed  plastic  tamgential  modulus  can  be  obtained 
by  eliminating  P*  from  Eqs  (45)  and  (67).  The  result 
gives 


=  2/iere, 


FJr  ~  (^1^1  +  4>2^2  +  4>i 

-b  <^2^2  +  £’’3^3 


•  (74) 


At  this  juncture  ,  em  additional  comment  about  P  is  nec- 
essaxy.  It  is  cleair  from  Eq  (30)  that  for  a  purely  viscoelas¬ 
tic  deformation  P  tends  to  zero.  If  p  becomes  00,  on  the 
other  hand,  the  deformation  has  no  apparent  hardening, 
i.e.,  becomes  zero  in  Eq  (45).  There  is  a  possibility 
that  shear  banding  may  occur  at  a  softening  stage,  i.e., 
when  ^  <  0.  In  other  words,  shear  banding  may  occur  at 
any  value  of  p.  Thus,  to  predict  the  onset  of  sheen  band¬ 
ing,  it  is  necessary  to  include  the  evolution  of  the  loading 
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Table  1;  Effect  of  a  on  shear  band  angles  9i  and  mea¬ 
sured  form  the  flow  direction. 


Q 

Or 

O2 

-0.080 

18.08 

103.94 

-0.090 

18.32 

103.69 

-0.100 

18.57 

103.45 

-0.110 

18.81 

103.21 

-0.118 

19.00 

102.99 

-0.120 

19.05 

102.97 

-0.130 

19.29 

102.73 

-0.140 

19.52 

102.50 

-0.150 

19.75 

102.27 

-0.160 

19.98 

102.04 

path  as  a  part  of  the  analysis  because  plastic  deformation 
is  path  dependent,  and  the  analysis  would  require  numer¬ 
ical  simulation.  This  is  beyond  the  scope  of  the  current 
paper,  and  will  be  addressed  in  a  later  paper. 

Orientation  of  shear  bands 

The  characteristic  directions  of  shear  bands,  in  general, 
are  given  by  Eq  (71)  from  which  the  orientation  of  shear 
bands  6  with  respect  to  the  directions  of  principal  stresses 
can  be  found  by 

e  =  ±tan-^A.  (75) 

In  order  to  compare  the  analytic  predictions  with  the  ex¬ 
perimental  results  provided  by  Bair  et  al.  (1993),  an  ap¬ 
propriate  coordinate  transformation  must  be  applied  to 
obtain  the  shear  band  angles,  &i  and  62,  which  are  mea¬ 
sured  from  the  solid  surface,  i.e.,  the  flow  direction,  in 
their  high  pressure  flow  visulization  cell. 

It  was  mentioned  earlier  that  the  value  for  a  used  in  the 
Drucker-Prager  plastic  dissipation  potential  is  a  material 
property,  and  needs  to  be  determined  from  experiments. 
To  determine  a  without  knowing  the  precise  loading  his¬ 
tory,  we  assume  that  l//3  =  0  because  ^er  is  in  the  order 
of  (Lee,  1989).  We  also  note  that  for  most  engi¬ 

neering  materials,  if  not  for  all,  W  is  relative  small  when 
compared  with  unity,  thus  1  -|-  jV  »  1.  The  above  two 
assumptions  together  with  Eqs  (50),  (67),  (69),  and  (71) 
are  used  to  determine  the  orientation  of  shear  bands  for 
different  values  of  a.  Using  0.4  for  v  and  0.5  for  the 
transverse  stress  ratio  k  which  is  defined  as  <73/(0’!  *f  (Tj), 
we  tabulate  the  results  in  Table  1.  It  is  to  be  expected 
that  a  would  have  an  effect  on  the  orientation  of  shear 
bands.  However,  as  shown  in  Table  1  the  effect  is  mod¬ 
est;  for  the  range  of  a  which  has  physical  significance  and 
which  covers  a  class  of  polymers  (Bowden,  1973),  the  ef¬ 
fect  of  a  is  less  observable.  Note  that  a  is  negative  in 
the  sense  that  the  shear  component  of  ^  increases  as  the 
pressure  incresses  when  <j>  is  taken  as  the  loading  func¬ 
tion  and  the  flow  rule  is  associative.  This  is  consistent 


with  experimental  observations  that  shear  strength  of  a 
lubricant  increases  as  pressure  increases.  These  values 
of  a  lie  within  the  range  of  vsdues  predicted  by  Drucker 
and  Prager  (1952).  For  polyethyl  ether  5P4E,  Bair  et. 
al.  (1993)  reported  shear  band  angles  of  61  =  19”  and 
$2  =  103°.  As  shown  in  Table  1,  their  result  is  matched 
by  the  case  for  which  a  =  —0.118.  Therefore,  we  use  this 
value  for  all  subsequent  calculations.  It  should  be  noted, 
however,  that  more  experimental  data  may  be  required 


-3  -2  -1  0  1  2  3  4  5 

FIG.  2.(a) 


<72/<Ti 

FIG.  2-(b) 


FIG.  2  Orientation  of  shear  bands  aX  v  =  k  =  0.4. 


to  corroborate  the  value  of  a. 

Figure  2  shows  shear  band  angle  as  a  function  of  the 
ratio  of  two  in-plane  principal  stresses,  <72/<ri,  which  is  an 
invariant.  This  result  is  plotted  for  Poisson’s  ratio  »'=0.4 
and  for  the  transverse  stress  ratio  i  =  0.4.  When  the  ra¬ 
tio  <72/t7i  takes  on  values  such  that  it  is  in  between  0.761 
and  1.313  as  shown  by  the  dotted  area,  the  field  equa¬ 
tions  are  elliptic  and  hence  no  shear  bands  will  occur. 
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FIG.  3-(a)  FIG.  4-(a) 
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FIG.  3  Orientation  of  shear  bands  at  =  0.4.  FIG.  4  Orientation  of  shear  bands  at  t'  =  0.45. 


In  general,  shear  bands  are  unlikely  to  occur  when  the 
loading  state  is  dominated  by  its  hydrostatic  component. 
Outside  the  shaded  region,  i.e.,  either  <r2/<ri  <  0.761  or 
02! (Ti  >  1.313,  the  field  equations  become  hyperbolic  and 
thus  line  or  characteristic  solutions  exist  and  she2u:  band¬ 
ing  becomes  possible.  This  result  is  repeated  in  Fig.  3 
for  a  further  comparison  with  the  other  two  cases  which 
have  different  values  of  k. 

For  a  liquid  lubricant,  the  loeuiing  state  excludes  the 
possibility  of  hydrostatic  tension,  thus  we  present  in  the 
following  the  study  of  shear  banding  for  (T2/<r\  >  1  that 
corresponds  to  XmasIP  >  0,  where  Tmax  =  (o'i  -  (^7)/^  is 
the  maximum  shear  stress  and  is  non-negative.  Note  that 
XmaxtP  is  also  2m  invariant  and  is  non-negative.  Figures 
3  to  5  show  the  orientation  of  shear  bands  under  different 
loading  states  for  which  the  Poisson’s  ratio  u  are  0.4, 0.45, 
and  0.499,  respectively.  For  each  case,  the  results  are 
shown  for  different  values  of  k  which  is  bounded  by  v  and 
0.5  in  a  plane  straun  deformation.  The  results  are  shown 


as  a  function  of  tmax/P-  As  an  example,  by  comparing 
the  data  for  k  =  0.45  between  Fig.  3  and  Fig.  4,  we 
find  little  effect  of  v  on  the  changing  orientation  of  shear 
bands,  at  least  for  the  scale  shown.  By  examining  the 
results  at  a  specific  value  of  Tmax  IP-,  we  in  fact  know 
that  there  is  a  slight  increase  in  61  and  a  slight  decrease 
in  ^2  as  1/  is  increased.  For  Tmax/P  —  0.105,  at  which 
the  shear  band  angles  were  measured  by  Bair  and  Winer 
(1993),  and  with  v  =  0.4,  it  can  be  observed  from  Fig. 
3  that  shear  banding  is  not  possible  at  ifc  =  0.4,  whereas 
shear  bands  occur  at  9x  =  19®  and  62  =  103®  for  k  — 
0.5.  The  effect  of  increasing  transverse  stress  ratio  on  the 
orientation  of  shear  bands  follows  the  same  trend  as  that 
of  V,  however,  the  changes  in  orientation  are  faur  more 
pronounced.  This  observation  of  the  effects  of  Poisson’s 
ratio  and  the  tramsverse  stress  ratio  on  shear  banding  is 
consistent  with  that  reported  by  Lee  (1989)  in  the  context 
of  elastoplasticity.  The  magnitude  of  k,  in  general,  varies 
during  the  course  of  deformation.  As  the  deformation 
becomes  more  amd  more  incompressible  in  nature,  the 
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magnitude  of  k  tends  to  its  limit  at  0.5.  For  a  nearly 
incompressible  medium,  the  inclination  of  the  shear  bands 
appears  to  be  independent  of  the  pressure.  This  can  be 
seen  &om  Fig.  5,  and  the  result  is  consistent  with  the 
observation  made  by  Bair  et  al.  (1993).  FVom  Figs  3  to 
5,  in  addition,  it  can  be  seen  that  the  angle  of  shear  bands 
appears  to  reach  a  limiting  value  as  the  deformation  tends 
towards  an  increasingly  incompressible  behavior.  To  be 
more  specific,  this  angle  is  between  19**  and  19.25**  for  $i 
and  between  102.77”  and  103”  for  $2,  depending  on  the 
magnitude  of  the  Poisson’s  ratio  v. 


"^mas/p 

FIG.  5-(b) 

F1G.5  Orientation  of  shear  bands  at  =  0.499. 

Another  observation  is  that  with  an  increasing  trans¬ 
verse  stress  ratio,  i.e.,  as  the  deformation  becomes  in¬ 
creasingly  incompressible  in  nature,  the  regime  of  elliptic- 
ity  or  the  range  of  stress  ratio  for  which  the  field  equations 
are  elliptic  becomes  smaller.  Eventually  at  the  stress  ra¬ 
tio  of  unity,  a  hydrostatic  state  of  stress  results,  and  no 
shear  bands  can  be  observed. 


CONCLUDING  REMARKS 

We  give  in  this  paper  a  constitutive  model  for  vis¬ 
cous  materials  which  have  both  rate-dependent  and  rate- 
independent  deformation  behavior  under  a  static  condi¬ 
tion.  Based  on  this  model,  the  analysis  also  focuses  on 
identifying  conditions  which,  if  achieved  in  a  deforma¬ 
tion,  would  produce  a  change  of  type  in  the  controlling 
field  equations.  The  underlying  hypothesis  is  that  shear 
bands  will  appear  when  the  character  of  the  field  equa¬ 
tions  changes  from  elliptic  to  hyperbolic.  The  results  are 
qualitatively  in  accord  with  experimental  observations. 
This  approach  to  shear  banding  problem  is  by  now  stan¬ 
dard  and  is  applicable  to  both  solid  and  liquid  media, 
at  least  within  the  framework  of  the  particular  decom¬ 
position  hypothesis  proposed  in  this  study.  It  should  be 
noted  that  there  is  by  no  means  imanimity  on  the  proper 
decomposition  of  viscoelastic  and  rate-independent  strain 
or  of  the  proper  relation  between  stress  and  strain  rates. 
Different  relation  would  give  results  that  may  be  differ¬ 
ent  in  some  respects.  In  addition,  since  the  evolution  of 
the  stress  path  is  unknown  in  advance,  one  should  note 
that  the  figures  in  the  last  section  can  only  be  used  qual¬ 
itatively  without  solving  an  actual  deformation  problem. 
Thus  one  should  also  note  that  all  parts  of  the  diagrams 
may  not  be  actually  accessible.  It  is  not  known  whether 
the  increased  and  decreased  angles  for  bands  within  a  very 
narrow  window  of  stress  ratios  would  be  predicted  for  any 
realizable  deformation  history.  Nevertheless,  the  experi¬ 
mental  data  on  the  fully  developed  shear  bands  does  give 
credence  to  the  theory  in  a  qualitative  sense.  Thus  the 
next  focus  of  research  should  be  a  numerical  implemen¬ 
tation  of  the  model  in  order  to  bring  out  the  full  features 
of  the  theory. 
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